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ABSTRACT 

The  cavitation  performance  of  a  seven-biaded,  highly 
skewed  nropeller  for  a  Naval  auxiliary  (AO-177)  is 
evaluated  using  lifting  surface  numerical  hydrodynamic 
methods.  An  important  contribution  to  accurate  predition 
of  cavitation  performance  is  shown  to  be  an  accurate  model 
of  the  effective  wa  k e . 

A  new  model  of  the  effective  waice  is  described.  A 
program  to  calculate  unsteady  time-averaged,  but  not 
circumferentialiy  averaged,  field  point  velocities  is 
presented.  The  propel ler-induced  velocities,  along  with 
the  original  nominal  wake,  are  combined  with  Huang's 
axisymmetric  effective  wake  scheme  in  pie-shaped  wake 
segments  to  determine  the  effective  wake. 

The  method  is  used  to  predict  severe  cavitation 
eatent.   This  prediction  is  confirmed  by  S3PA  experiments. 
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INTRODUCTION 


The  David  W.  Taylor  Naval  Ship  Research  and 
Development  Center  has  recently  designed  a 
seven-b 1 aded  ,  highly  skewed  propeller  for  a  Naval 
Auxiliary  (AO-177).  This  design  is  instructive 
concerning  the  state  of  the  art  in  computer-aided 
propeller  design. 

The  history  of  the  propeller  design  is  chronicled 
in  reports  by  Boswell  (1),  Valentine  and  Chase  (2),  and 
Hendrican  and  Remmers  <3) .  Computer-aided  propeller 
applications  included  lifting  line  theory,  used  to 
determine  the  radial  load  distribution  and  the  radial 
hydrodynamic  pitch  angle,  and  lifting  surface  theory, 
used  to  determine  the  final  geometric  pitch 
distribution,  camber  distribution,  and  f inai  propeller 
offsets,  including  fillets,  trailing  and  leading  edge 
details,  additional  thickness  added  to  the  trailing 
edge,  tip  geometry,  rake,  and  hub  details.  Cavitation 
performance  was  predicted  by  the  method  of  Burr  i  I  I  and 
Emerson  (4)  and  was  also  evaluated  by  experiment. 

However,  sea  trials  of  the  AO-177  indicated  a 
severe  unsteady  propeller  force  problem.  Inspection  of 
the  propeller  indicated  that  it  was  cavitating 
significantly  worse  than  had  been  predicted.  The 
Burr  i  I i   and  Emerson   method,    originally   designed   for 
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f our-t o-si  n-b laded  propellers  with  an  expanded  area 
ratio  of  0.60,  did  not  extend  well  to  the  seven-b laded 
AO-177  propeller  with  an  expanded  area  ratio  of  0.77, 
even  though  both  were  for  merchant-type  hulls. 
Further,  scale  effects  could  have  thrown  off  the  model 
tests  . 


/ 
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II.    IMPACT  OF  EFFECTIVE  WAKE 


The  AO-177  propeller  performance  was  evaluated  at 
MIT  using  some  numerical  propeller  analysis  computer 
programs  developed  at  MIT.  Kerwin  and  Lee  (5)  describe 
MiT-PUF-2,  a  computer  program  to  calculate  forces, 
steady  and  unsteady,  generated  by  a  propeller. 
M1T-FPV,  a  steady  field  point  velocity  program,  is 
described  in  a  report  by  Min  (6).  Lee  (7)  reports  on 
M1T-PUF-3,  a  computer  program  to  predict  steady  and 
unsteady  propeller  cavitation.  extent  and  volume.  All 
three  programs  .  model  the  propeller  as  a  grid  of 
discrete  vortex  segments.  The  wake  is  also  modelled  as 
a  vortex  grid  in  each  program. 

Cavitation  performance  calculations,  more  than 
propeller  force  and  induced  velocity  calculations,  are 
highly  sensitive  to  the  wake  data  used  as  input. 
Experience  at  MIT  and  elsewhere  has  shown  that  using 
the  nominal  wake  data  to  approximate  inflow  conditions, 
wh  ile  acceptable  for  propeller  force  and  induced 
velocity  calculations,  is  not  acceptable  for  cavitation 
performance  calculations.  The  wake  velocities  measured 
behind  a  model  in  the  absence  of  any  propeller  (the 
nominal  wake)  must  be  modified  to  account  for  change  in 
boundary  conditions  imposed  by  the  presence  of  a 
propeller.  This  new,  modified  wake  is  generally  called 
the   "effective   wake."    There   is,   unfortunately,   no 
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general  method  for  anticipating  the  effective  wake 
generated  by  a  particular  propeller  in  a  general 
nomi  na  1  wa  Ic  e  . 

Huang  and  others  <8,9>  have  developed  both  theory 
and  numerical  schemes  for  calculating  the  effective 
wake  of  a  propeller  in  an  asisymmetric  nominal  wake. 
Det  Norske  Veritas  has  attemoted  to  adapt  Huang's 
method  to  non-as 1 s ymme t r i c  wakes,  such  as  surface  ship 
wakes,  by  dividing  the  nominal  wake  into  pie-shaped 
segments  and  then  applying  the  asisymmetric  wake 
calculation  within  each  pie  segment,  using  the  nominal 
wake  and  the  steady  propeller  induced  velocities.  This 
approach  is  an  improvement  over  blindly  applying  the 
asisymmetric  effective  wake  calculation  to  a  surface 
ship's  wake,  but,  owing  to  the  very  sharp  AO-177  wake 
(even  by  surface  ship  standards),  it  was  judged  that 
some  further  refinement  of  Huang's  method  was 
r  equ  i  r  ed . 

Huang's  method  assumed  an  asisymmetric  nominal 
wake  and  a  steady  propeller-induced  velocity  field. 
The  Det  Norske  Veritas  method  assumed  a  general  nominal 
wake  and  a  steady  propel ler-induced  velocity  field. 
Short  of  a  completely  general  effective  wake 
calculation,  the  most  obvious  refinement  was  to  assume 
a  general  nominal  wake,  but  allow  an  unsteady 
propel  ler-induced   velocity  field.   This  was   the  method 
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used  for  this  paper  . 

The  current  method  of  effective  wake  calculations 
can  be  summarized  in  the  following  steps: 

1)  Use  the  nominal  wake  and  propeller  data  to 
calculate  the  unsteady  propeller-induced  velocity  field 
at  representative  points  both  inside  and  outside  the 
propeller  disk. 

2)  Divide  the  nominal  wake  into  small  pie-shaped 
segments  as  in  the  Det  Norske  Veritas  method. 

3)  Calculate  the  effective  wake  velocities  within 
each  pie  segment  using  the  nominal  wake  velocities  and 
the  unsteady  p r o p e  1  1  e r - i nduc e d  velocities,  assuming 
that  each  pie  segment  acts  as  though  it  were  part  of  an 
asisymmetric  wake  with  the  same  radial  velocity 
distribution. 


■y 
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III.    METHOD  OF  EFFECTIVE  WAKE  CALCULATION 


The  two  new  aspects  of  the  current  effective  wake 
caiculations  are  the  propeller-induced  unsteady  field  point 
velocity  (UFPV)  calculation  and  the  use  of  these 
time-averaged  out  not  c i r c umf e r en t  i a  1  1 y -a v e r a g e d  velocities 
in  9  pie  segment  nominal  wake  modification  calculation 
( P I EWAKE ).  Each  of  these  two  new  methods  will  be 
illustrated  by  a  typical  calculation. 

The  unsteady  field  point  velocity  (UFPV)  program  is 
basically  a  generalization  of  the  steady  field  point 
velocity  program  (MIT-FPV)  by  Min  .  <6>.  The  unsteady 
version  uses  the  full  wake  lattice  arrangement  also 
employed  in  MIT-PUF-2.  Program  inputs  are  the  FILE14  DATA 
file  output  from  MIT-PUF-2  and  interactive  directions. 
Program  outputs  are  either  single  point  calculations 
directed  to  the  user's  terminal  or  a  data  file  to  be  used 
as  input  into  PIEWAKE,  for  effective  wake  caiculations. 

UFPV  first  "reads  the  FILE  14  DATA  input,  containing 
mostly  propeller  geometry,  singularity  strengths,  and 
singularity  geometry.  It  should  be  noted  that  the  number 
of  propeller  revolution  time  steps  specified  in  MIT-PUF-2 
must  be  divisible  by  the  number  of  propeller  blades  For 
effective  wake  caiculations,  the  wake  input  to  MIT-PUF-2 
ought  strictly  to  be  the  effective  wake;  however,  using  the 
nominal  wake  to   calculate  the   propel  ler-induced  velocities 
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is  a  good  approximation. 

The  user  next    inputs   field   point  coordinates.     The 

field    point  position   may   be   entered     in    cylindrical 

coordinates  or  in  terms  of    the  rotating  propeller   and  wake 
grids  . 

t 

The  velocity  induced  by  one  blade  at  each  time  step  is 
then  calculated.  At  each  time  step,  the  vortex  lattice  is 
assigned  the  strengths  corresponding  to  the  time  step  of 
the  propeller  rotation.  The  velocity  induced  by  each 
element  of  the  vortex  grid,  both  blade  and  wake,  is  tnen 
summed  to  obtain  the  total  velocity  induced  by  one  blade  in 
that  particular  angular  position. 

The  velocities  induced  by  each  blade  of  the  propelier 
are  then  summed  to  produce  the  total  propeller-induced 
velocity  for  that  position  of  the  propeller.  Velocities 
are  calculated  for  different  propeller  positions  and  then 
harmonically  analyzed.  The  zeroeth  harmonic,  then,  is  the 
time-averaged  propeller-induced  velocity  for  that  field 
point.  <By  averaging  all  points  at  a  given  radius  we  could 
obtain  the  c i r c umf e r en t  i a  1  1  y - a v e r a g e d  ,  time-averaged, 
propeller-induced  velocity,  the  output  of  the  steady 
M1T-FPV . ) 

For  effective  wake  calculations,  such  velocities  are 
computed  for  up  to  60  points  per  radius  for  radii  varying 
from   propeller  hub   to   1.7   times   the   propeller   radius. 
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Between  the  propeller  hub  and  propeller  tip,  the  velocities 
are  calculated  at  the  leading  edge  panel  on  the  propeller 
grid. 

For  the  AO-177  propeller  analysis,  56  points  per 
radius  for  eight  radii  were  used  to  create  a  file  of 
unsteady  field  point  velocities. 

PlEVAKE  is  a  program  to  perform  Huang's  axisymmetric 
effective  wake  contraction  calculations  in  each  of  many 
pie-shaped  segments  of  the  wake.  It  uses  nominal  wake  data 
and  unsteady  propeller-induced  velocity  da: a  from  U F  P  V  as 
inputs  to  produce  an  effective  wake  velocity  data  file  as 
output  . 

After   reading   nominal    wake   data   and  the'   induced 

velocity   field  data,  PIEVAKE  extrapolates  the  nominal  wake 

data  to   the  hub.     The   innermost   radius   of  the   induced 

velocity   field   is  taken  as   the  hub.    (UFPV  had  taken,  as 

its  innermost  radius,  the  hub  radius  used  in  MIT-PUF-2  and 
passed  to  UFPV  by  way  of  FXLE14  data  file.) 

In  each  pie  segment,  Huang's  method  is  used  to 
calculate  an  effective  velocity  and  the  affective  radius 
corresponding  to  that  effective  velocity.  In  the  numerical 
form,  this  method  is  iterative  in  nature. 

A  finite  difference  equation,  given  in  reference  (9), 
is  used  to   calculate  the   effective   velocity.     Given   the 
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nominal  wake  at  radii  (assumed,  for  the  first  iteration,  to 
be  the  nominal  radii)  in  the  presence  of  a  propeller,  an 
effective  velocity  can  be  calculated  if  an  effective 
velocity  at  the  next  radius  away  from  the  hub  is  known. 
The  effective  velocity  at  the  outermost  radius  is  assumed 
to  be  equal  to  the  nominal  velocity  there,  enabling  all 
effective  velocities  to  be  calculated. 

A  second  finite  difference  equation,  also  given  in 
reference  (9),  is  used  to  calculate  the  effective  radius 
corresponding  to  the  set  of  effective  velocities  just 
calculated.  This  equation  uses  nominal  radii,  nominal 
velocities,  and  apparent  velocities.  (An  apparent  velocity 
is  the  sum  of  the  effective  velocity  and  the 
propeller-induced  velocity  at  a  point.)  Given  one 
effective  radius,  the  next  radius  away  from  the  hub  can  be 
calculated.  Since  the  innermost  radius  of  all  velocity 
fields  is  set  at  the  hub  radius,  all  effective  radii  can  be 
determined  sequentially  from  innermost  to  outermost.  These 
radii  can  then  be  used  in  the  finite  difference  equation  to 
calculate  the  neat  iteration  of  effective  velocities. 
PIEWAKE  iterates  five  times  for  each  pie  segment. 

In  addition  to  the  effective  velocity  wake  field, 
PIEVAKE  also  computes  the  mean  velocity,  the  volumetric 
mean  velocity  of  the  nominal  wake,  the  volumetric  mean 
velocity  of  the  effective  wake,  and  the  effective  blockage. 

The  effective   velocity   field  then   then   be  processed 
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and  used  by  MlT-PUF-3  as  any  other  waice  field. 
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IV.    AO-177  EFFECTIVE  WAKE  CALCULATIONS 


Appendix  I  contains  plots  of  nominal,  induced,  and 
effective  velocities  varying  circumferential I7  at  several 
specified  radii.  Appendix  II  contains  the  same 
information,  except  that  the  nominal  and  effective 
velocities  are  normalized  on  the  respective  volumetric  mean 
velocity.  These  plots  confirm  that  the  modifications  to 
the  nominal  wake  outlined  in  Chapter  III  are  indeed 
reasonable.  Several  salient  conclusions  can  be  drawn  from 
these  plots: 

1)  Effective  velocity  modifications  of  the  nominal 
velocities  are  in  the  expected  direction.  Physically,  the 
imposition  of  the   propeller  on  the  nominal   wake  draws  more 

water  into  the  slipstream,  thus  increasing  the  velocities. 
This  "faster  water"  will  come  from  radii  greater  than  the 
field  point  being  sampled;  in  other  words,  the  wake  will 
seem  to  "contract." 

2)  Effective  velocity  modifications  of  the  nominal 
wake  velocities  are  greatest  at  lower  radii.  The  wake 
"contraction"  will  seem  to  concentrate  here. 

3)  In  general,  the  wake  peak  is  more  narrow  in  the 
effective  wake  than  in  the  nominal  wake.  In  the  case  of 
the  AO-177  nominal  wake,  the  existence  in  the  presence  of  a 
propeller  of  a  wake  peak  as  sharp  as  that  in  the  nominal 
wake  is  counter-intuitive.  Therefore,  a  decrease  in  the 
wake  peak  is  to  be  expected. 

4)  The  normalized   effective   velocities  at  0   degrees 
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are  not  appreciably  different  from  that  of  the  nominal 
wake.  This  is  disappointing  and  was  not  expected.  It  may 
be  a  result  of  radial  vorticity,  ignored  by  the 
axisymmetric  effective  wake  calculation  performed  in  each 
pie  segment.  This  assumption  is  least  appropriate  in  the  0 
degree  region  of  the  wake.  A  more  general  effective  wake 
calculation  scheme  may  be  required  to  obtain  better  results 
in  this  region. 

Results  of  M1T-PUF-3  cavitation  performance  of  the 
AO-177  propeller  are  shown  in  Appendices  III  and  IV.  As  an 
indicator  of  the  influence  of  effective  wake  modifications 
to  the  nominal  wake,  MIT-PUF-3  was  run  twice,  once  with  the 
nominal  wake  as  input  and  once  with  the  effective  wake  (the 
more  accurate  boundary  condition)  as  input.  Appendix  III 
shows  cavitation  performance  given  the  nominal  wake  as 
input  while  Appendix  IV  shows  cavitation  performance  given 
the  effective  wake  as  input.  In  general,  the  plots  show 
that  the  cavitation  extent  is  not  changed  much,  although 
cavitation  performance  is  better  in  the  effective  wake  at 
the  lower  radii  where,  as  has  been  shown,  the  effective 
wake  modifications  are  the  greatest.  At  the  outer  radii, 
where  the  effective  wake  modification  is  the  least 
significant,  the  cavitation  performance  is  least  changed 
between  the  nominal  wake  and  effective  wake  cases. 

Appendix  V  contains  a  graph  comparing  cavitation 
extent  of  MIT-PUF-3  output  given  nominal  wake  input  against 
that   given   effective   wake  input.    It  plots  the   length  of 
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cavitation  on  the  seventh  radial  cavitation  panel  (about 
the  .8  propeller  radius)  non-d i mens i ona 1 i z e d  on  the  local 
chord  length  as  the  blade  rotates  through  360  degrees. 
This  graph  indicates  that  the  effective  wake  modification 
did  not  reduce  cavitation  extent  on  a  blade,  a  result  noted 
earlier.  However,  a  larger  contrast  between  the  nominal 
waJce  and  effective  wake  cases  is  shown  in  Appendix  VI  which 
graphs  the  cavitation  volume  on  a  blade  in  the  two  cases. 
The  volume  is  non-dimensionalized  on  the  cube  of  the 
propeller  radius.  This  graph  shows  a  drop  of  approximately 
35 'A  of  peak  cavitation  volume.  This  change  in  cavity 
volume  demonstrates  the  sensitivity  of  cavitation 
performance  to  small  changes  in  the  wake  profile.  This 
confirms  that  efforts  to  refine  effective  wake  calculations 
were  indeed  justified. 

The  effective  wake  cavitation  performance  indicated  by 
M1T-PUF-3  compares  reasonably  well  with  experiments 
performed  by  the  SSPA  propeller  tunnel.  These  experiments 
measured  cavitation  performance  of  the  A  0  -  1  7  7  propeller  by 
photographing  model  tests.  While  no  cavitation  volume 
results  were  obtained,  cavitation  extent  could  be  compared 
with  M1T-PUF-3  oredictions.  Despite  reductions  of  the 
effective  wake  calculations  compared  with  nominal  wake 
calculations,  MIT-PUF-3  still  tended  to  predict  greater 
cavitation  than  the  experiments  indicated,  particularly  at 
the  middle  radii.  This  seems  to  indicate  that  the 
effective  wake  modifications,  while  tending  to  the  proper 
direction,  were  not  large   enough.   This  is  consistent   with 
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the  conclusions  from  comparing  the  two  wakes  that  a  more 
general  effective  wa k e  scheme  including  radial  vorticity  in 
high  wake  regions  is  required. 


-  18  - 


CONCLUS IONS 


1)  Calculation  of  an  effective  wake  using  Huang's 
axisymmetric  method  in  pie-shaped  segment  with  unsteady 
propel Isr-induced  velocities  is  a  reasonable  approximation. 
Tt     still    tends   to    underestimate   the   effective    wake 

modification  to  the  nominal  wake  in  the  vicinity  of  0 
degrees.  To  correct  this,  a  more  general  effective  wake 
calculation  scheme  including  radial  vorticity  effects  may 
be  required. 

2)  The  current  effective  wake  calculation  scheme  can 
be  used  to  predict  cavitation  performance  for  highly  skewed 
propellers  with  up  to  seven  blades. 

3)  Cavitation  performance  is  indeed  sensitive  to 
small  changes  in  the  wake  profile.  Seemingly  small 
improvements  in  effective  wake  calculations  yield 
significantly  more  accurate  cavitation  performance 
predict  ions  . 
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VI.    RECOMMENDATIONS 


1)  A  more  general  effective  wake  calculation  scheme, 
including  the  effects  of  radial  vorticity,   is  required. 

2)  The  computational  time  (and  hence  cost)  of  the 
unsteady  field  point  velocity  program  (UFPV)  can  probably 
be  cut  without  significant  sacrifice  of  accuracy  by 
decreasing  the  number  of  points  calculated  per  radius. 
Appendix  I  shows  that  the  variation  of  induced  velocity 
with  rotation  is  slight  and  smooth.  Calculation  of  fewer 
points  and  estimation  of  the  remaining  points  using 
harmonic  analysis  would  probably  not  reduce  accuracy. 
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**  APPENDIX  I  ** 


Nomina.!,   Induced,  and  Effective  Velocities 
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**  APPENDIX  II   ** 

Nominal,   Induced,  and  Effective  Velocities 
Non-Dimens  ional  i  zed 
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**  APPENDI X  III   ** 


Propeller  Cavitation,  Nominal  Va k e 
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**  APPENDIX  IV  ** 


Propeller  Cavitation,  Effective  Wak 


/ 


44  - 


© 

o 

id  o 

> 

+» 

o 


-  45  - 


o 

3: 
O 

> 


o 

«m 


O 


C\! 


CM 


ON 


-  46  - 


8 

:* 
o 

©   CM 

>     T-l 

•H 

-P- 
O 
6 


ta 


-  hi  - 


a 

o 

O  CO 

>   tH 


o 


-  48  - 


> 

-P 
O 

CD 
Cm 

u 


CM 


-  k9  - 


0) 

> 

-P 
O 

<u 


o 
en 


o 

co 
CM 


en 


-  50  - 


o 

o 
>  r^ 


o 


-  51  - 


o 

0)  CM 
>  J" 


O 
CD 

Cm 


-  52  - 


**  APPENDIX  V  ** 


Cavity  Length,  Nominal  and  Effective  Wakes 
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**  APENDIX  VI   ** 


Cavity  Volume,  Nominal  and  Effective  Wakes 
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**  APPEND  I  X  VII   ** 


SSPA  Exner  imental  Prooelier  Cavitation  Extent  Plots 
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**  APPEND  IX  VIII   ** 


Nominal  Walce  Diagram 
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**  APPENDIX  IX 


Etfective  Wake  Diagram 
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**  APPENDIX  X  ** 


UFPV  Program  Listing 
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c 
c 
c 

c 
c 
c 

c 

c 
c 

c 
c 
c 
c 

c 

c  »« 

c 

905 


MAIN  PROGRAM  IN  FPUV.  A  PROGRAM  TO  CALCULATE 
UNSTEADY  VELOCITY  HARMONICS  AT  FIELD  POINTS 
SPECIFIED  INTERACTIVELY  IN  RDTERM. 


»•■■»«»»»**•»»»*«»*»«»««•»««*: 


»»■**■*»***•»«  ..».«• 


*-«  COMMON  TO  CONTAIN  GEOMETRIC  PARAMETERS 


30 


COMMON  /  GEOCOM  /  MM , MN , NN , NW , NTORL . X( 1 2 , 10 ) . Y ( 1 2 . 10 ) . Z ( 1 2 . 10) , 

5  XTIP(  10.  1  1  ). YTIP(  10.  1  1  )  ,ZTIP(  10,  1 1 ) ,XH(  101 .7) ,YH(  101 ,7) . 

6  ZH(  101  .7) ,XW(2  1 .  10) . YW(21 .  10) ,ZW(21 .  10) . RH 

■«  COMMON  TO  CONTAIN  SINGULARITY  STRENGTHS 

COMMON  /  SINCOM  /  GT ( 10 . 10. 60) . G8 ( 10 . 9 . 60) . GTW( 10) . SB ( 10 . 9 ) , GTV 
&  . ,GAM(9.60) .GMEAN(9) . DKFW( 20 ) . HBWX . HBWY . HBWZ 

■»  COMMON  TO  CONTAIN  FIELD  POINT  DATA.  NBLADE.  IDENT.  NSR ,  AND  VSR 

COMMON  /  FPCOM  /  NBLADE . NSR . XFP . RFP . TZFP . YBAR . ZBAR . 
&    DBLADE,DTFP.-IDENT(  18)  ,TFP(60)  .VSR.  I  PLOT,  D I  AM 

■*  COMMON  TO  CONTAIN  VELOCITY  DATA 

COMMON  /  VELCOM  /  VIX ( 60 ) . VIT ( 60 ) . VIR ( 60) , U( 60 . 3 ) . NSAMP 
<*  COMMON  TO  CONTAIN  HARMONICS  COEFFICIENTS 

COMMON  /  HARCOM  /  A (  1 5 . 3  )  , 8 (  1 5 , 3 ) , AMP (  1 5 . 3 ) . PH(  1 5 . 3 ) . NH 
-*  ISTOP  IS  A  FLAG  SET  INTERACTIVELY  IN  ROTERM  TO  INOICATE  END  OF  P 

IST0P=O 

'*  RDFILE  WILL  READ  DISK  DATA  FILE  CREATED  BY  PUF-2.   NOTE  THAT  NSR 
PUF-2  INPUT  MUST  BE  DIVISIBLE  8Y  THE  NUMBER  OF  BLADES.  N8LADE. 

CALL  RDFILE 

-  INTERACTIVELY  READ  PARAMETERS  FOR  TERMINAL  SESSION 

FORMAT(/'  ENTER  NBLADE,  NTORL.  IPLOT.'/ 
&   '  FOR  HELP.  ENTER  ZERO  FOR  N8LADE.') 
WRITE(6,905) 

REAO( 5, •)  NBLAOE, NTORL. IPLOT 
IF  (NBLADE. LT.O)  GOTO  20 
IF  (NBLADE. EO.O)  GOTO  90 
IF  (IPLOT. EO. 2)  CALL  P I EDAT( I  STOP ) 
IF  (ISTOP. EO. 1)  GOTO  20 


UFP0OO10 
UFPOOO20 
UFP0O030 
UFP0O040 
UFPOOO50 
UFP0OO60 
UFP0O070 
UFPOOO80 
UFP0O090 
UFP00100 
UFP001 10 
UFP00120 
UFP00130 
UFPOO140 
UFP00150 
UFP0O160 
UFP00170 
UFP0O180 
UFP00190 
UFP00200 
UFP0O210 
UFP00220 
UFPOO230 
UFP00240 
UFP0O250 
UFP00260 
UFP00270 
UFP00280 
UFP00290 
UFP00300 
UFP0O310 
UFP00320 
UFPOO330 
R0UFPOO34O 
UFP0O350 
UFP0O360 
UFPCO370 
IUFP00380 
UFP00390 
UFP0O40O 
UFP004 10 
UFP00420 
UFP00430 
UFPOO440 
UFP0O450 
UFPOO460 
UFP00470 
UFPOO480 
UFP00490 
UFPOO50O 
UFP0O5  10 
UFP0O520 
UFP0O530 
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C  »**  RDTERM  GETS  DATA  INTERACTIVELY  FROM  TERMINAL  ANO  COMPUTES 

C      NECESSARY  FIELD  POINT  GEOMETRY. 

C 

10   CALL  RDTERM( ISTOP) 
C 

C  **«  IF  OPERATOR  INDICATES  END  OF  PROGRAM.  ISTOP  WILL  BE  SET  TO  1 
C 

IF  (ISTOP. GT.O)  GOTO  20 
C 

C  •**  RUN  THROUGH  FPSTEP  FOR  EACH  "TIME  STEP.' 
C 

DO  100  K=1 .NSR 
C 

C  ...  FPSTEP  CALCULATED  VELOCITY  INDUCED  AT  FIELD  POINT  BY  ONE  BLADE 
C      AT  "TIME  STEP"  K 
C 

10O   CALL  FPSTEP(K) 
C 

C  **•  SUMVEL  WILL  SUM  FIELD  POINT  VELOCITIES  INOUCED  BY  EACH  BLADE. 
C 

CALL  SUMVEL 
C 

C  -»*  FOUCAL  WILL  GENERATE  HARMONIC  COEFFICIENTS  FOR  INOUCED  FIELD  POI 
C      VELOCITY  AS  PROPELLER  ROTATES. 
C 

CALL  FOUCAL 
C 

C  -«-  PLOTVL  PLOTS  INDUCED  FIELD  POINT  VELOCITIES  AS  A  FUNCTION  OF  BLA 
C      ROTATION. 
C 
901   FORMAT('  AXIAL  COORDINATE '. 4X . F5 . 3/ '  RADIAL  COORD INATE ' , 2X . F6 . 3/ 
&   '  ANGULAR  COORDINATE'  .2X.F5.  1///'  AXIAL  HARMONI C '  . 6X , F8 . 5/ 
&   '  RADIAL  HARMONIC .5X.F8. 5/'  TANGENTIAL  HARMONIC  '.F8.5) 
WRITE (6. 901)  XFP.RFP.TZFP, (AMP( 1 ,U) . J=1 .3) 

IF  (IPL0T.EQ.1)  CALL  PLOTVL ( N8LADE . NSR ,NH . TFP , TZFP . U . A , B . AMP . PH, 
&   XFP.RFP.IDENT) 
GOTO  10 
906   FORMAT(/'  NBLAOE  <  0  =>  STOP' 

&        /8X.      '-  0  =>  HELP'/8X.'>  0  =>  NUMBER  OF  BLADES'/ 
&   '  NTORL   =  0  =>  TOTAL  EF F ECT ' /8X , ' =  1  =>  THICKNESS  ONLY'/ 
&   8X.'=  2  =>  LOADING  ONLY'/'  IPLOT   =  0  =>  DATA  TO  TERMINAL'/ 
&   8X.'»  1  =>  PLOT  --  NOT  IMPLEMENTED  05/02/82  --  RLJ' 
&   /8X.'=  2  =>  CREATE  PIEWAKE  DATA  FILE'/) 
90   WRITE (6. 906) 

GOTO  30 
20   STOP 
END 


UFP00S40 
UFP0O550 
UFP00560 
UFP0O570 
UFP00580 
UFP0O590 
UFP0O6CO 
UFPOO6  10 
UFP00620 
UFPOO630 
UFPOO640 
UFP00650 
UFP0O660 
UFP00570 
UFP0O680 
UFPOO690 
UFP007CO 
UFP0O710 
UFP00720 
UFPOO730 
UFP0O740 
UFP0O750 

NTUFP00760 
UFP0O770 
UFPOO780 
UFP00790 
UFP0O800 

DEUFPOO810 
UFP0O820 
UFP0O830 
UFP00840 
UFP00850 
UFP00860 
UFP0O870 
UFP0O880 
UFP00890 
UFP0O90O 
UFP00910 
UFP00920 
UFP00930 
UFP00940 
UFP0O950 
UFP00960 
UFP00970 
UFP00980 
UFP0O990 
UFP01OOO 
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SUBROUTINE  FPVEL ( K , XF . YF , ZF . VIX , VI Y . VIZ )  UFP01010 

C  UFP01020 

C  •*»  COMMON  TO  CONTAIN  GEOMETRIC  PARAMETERS  UFP01030 

C  UFP01040 

COMMON  /  GEOCOM  /  MM . MN . NN . NW , NTORL . X ( 1 2 . 10 ) , Y ( 1 2 . 10 ) , Z ( 1 2 . 10 ) ,    UFP01050 

&    XTIP(  10.  1 1 ) . YTIP( 10. 1 1 ) ,ZTIP( 10.  1 1 ) .XH( 101 .7) ,YH(  101 .7) ,  UFP01060 

&    ZH( 101 .7) ,XW(21 . 10) . YW(21 . 10) .ZW(21 . 10) ,RH  UFP01070 

C  UFP01080 

C    •*•    COMMON    TO    CONTAIN    SINGULARITY    STRENGTHS             ■     •  UFP01090 

C  UFP01 100 

COMMON  /  SINCOM  /  GT( 10. 10. 60) . GB ( 10. 9 . 60) . GTW( 10) . S3 ( 10 . 9 ) . GTV    UFP01110 

&  ,GAM(9.60) ,GMEAN(9) .DKFW(20) . HBWX .HBWY . HBWZ               •  UFP01120 

C  UFP01130 

C  •**  COMMON  TO  CONTAIN  FIELD  POINT  DATA.  NBLADE.  IDENT.  NSR.  AND  VSR    UFP01140 

C             '  UFP01  150 

COMMON  /  FPCOM  /  NBLADE . NSR . XFP . RFP . TZFP . YBAR . ZBAR .  UFP01160 

&    DBLAOE.DTFP, IDENT(  18)  ,TFP(60) .VSR, IPLOT.DIAM  UFP01170 

C  UFP01  180 

DIMENSION  UTX(20.2).UTY(20.2) ,UTZ(20.2) . UBX ( 20) . UBY ( 20) .UB2(20)    UFP01 190 

&  ,VOFWX(20.9).VOFWY(20.9) ,V0FVZ(20.9)  UFP01200 

C  UFP01210 

C  •••  INITIALIZE  VARIABLES  UFP01220 

C  UFP01230 

MID=MN/2  ' UFP01240 

NVFW=NW  UFP01250 

UMPFW=1  UFP01260 

WAKX=0.  UFP01270 

WAKY=0.  UFP01280 

WAKZ=0.  UFP01290 

VIX=0.0  UFPO130O 

VIY=0.0  UFP01310 

VIZ=0.0  UFP01320 

DO  1  M=1,MM  UFP01330 

DO  1  N=1.NN  UFP01340 

C CHORDWISE  VORTICES  ON  BLADE UFP01350 

IF(NT0RL.E0. 1 )  GO  TO  8  UFP01360 

CALL  VORSEG(XF.YF.ZF.X(N.M) ,Y(N.M).Z(N.M).X(N+1.M).Y(N+1,M).  UFP01370 

1   Z(N+1 ,M),CVX.CVY.CVZ.SX.SY,SZ.O)  UFP01380 

C SPANWISE  VORTICES  ANO  SOURCES  ON  BLADE - UFP01390 

8  CALL  VORSEG(XF.YF.ZF,X(N.M) . Y(N.M) ,Z(N.M) ,X(N.M+ 1 ) . Y(N.M+1 ) .  UFP01400 
1   Z(N.M+1) .SVX.SVY.SVZ.SX.SY.SZ. 1 )  UFP01410 

SBNM=S8(N.M)  UFP01420 

G3NM=GB(N.M,K)  UFP01430 

GTNM=GT(N.M.K)  UFP01440 

IF(NTORL.EO.O)  GO  TO  9  UFP01450 

IF(NTORL.EO. 1 )  G8NM=0.0  UFP01460 

IF(NTORL. EO. 1 )  GTNM=0.0  UFP01470 

IF(NT0RL.E0.2)  S8NM=0.0  UFP01480 

9  VIX=VIX+GBNM*SVX+SBNM»SX+GTNM-CVX  UFP01490 
VIY»VIY+GBNM*SVY+S8NM»SY+GTNM»CVY  UFP015CO 
VIZ'VIZ+GBNM*SVZ+SBNM-SZ+GTNM-CVZ  UFP015  10 

1     CONTINUE  UFP01520 

C SEPARATED  TIP  VORTEX  ON  BLADE-- - UFP01530 
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.XTIP(N.LL) , 
,L1).ZTIP(N, 


35 


11 


73 
72 


IF(NTORL.EO.  1  )  GO  TO  10 

00  2  N=1 ,NN 

CVX=0.0 

CVY=0.0 

CVZ=0.0 

NTEM=NN-N+1 

DO  3  LL=1 .NTEM 

L1=LL+1 

CALL  VORSEG(XF.YF.ZF 
•   ,XTIP(N,L1).YTIP(N 

CVX=CVX+VX 

CVY=CVY+VY 

CVZ=CVZ+VZ 

VIX=VIX+G3(N.MM,K)*CVX 

VIY=VIY+G8(N,MM,K)*CVY 

VIZ=VIZ+GB(N,MM,K)*CVZ 
-TRANSITION  WAKE 

WAKX=0. 

WAKY=0. 

WAKZ=0. 

DO  35  N=1 , 100 

CALL  VORSEG(XF 
1   ZH(N+1  .  1  ) .VX 

WAKX=WAKX+VX 

WAKY=WAKY+VY 

WAKZ=WAKZ+VZ 

CALL  VORSEG(XF 
&   HBWX.HBWY.HBWZ 

HBWX=-HBWX 

HBWY=-HBWY 

HBWZ=-HBWZ 

DO  7  1  M=1 .MM 

IF(M.EO. 1 )  GO  TO 

DO  73  N=1 ,NW 

UTX(N. 1 )=UTX(N.2) 

UTY(N. 1)=UTY(N.2) 

UTZ(N, 1)=UTZ(N.2) 

DO  74  N=1 ,NW 

CALL  VORSEG(XF. YF.ZF.XW(N.M) , YV(N.M) , ZW(N.M) 
1   ZW(N,M+1 ) ,UBX(N) .UBY(N) .UBZ(N) .SX.SY.SZ.O) 

IF(M.GT..1  )  GO  TO  75 

CALL  VORSEG(XF,YF,ZF.XW(N,M) , YW(N, M) . ZW( N , M) 
1   ZW(N+1 .M) ,UTX(N, 1) .UTY(N. 1 ),UTZ(N. 1 ) . SX.SY 


YTIP(N.LL) .ZTIP(N.LL) 

L1 ) .VX.VY.VZ. SX.SY.SZ.O) 


YF.ZF.XH(N, 1 ) ,YH(N. 1 ) ,ZH(N, 1 ) 
VY.VZ. SX.SY.SZ.O) 


XH(N+1 .  1  )  . YH(N+1 .  1  ) 


YF,ZF.XH( 1 . 1 ),0. 
SX.SY.SZ.O) 


72 


.0. ,XH( 101. 1 ) .0. .0. 


XW(N.M+1 ) . YW(N.M+1 ) . 


XW(N+1 

sz.o) 


M)  ,YW(N+1 ,M) . 


75  CALL  VORSEG(XF. YF . ZF , XW(N. M+ 1 ) . YW( N. M+ 1 ) . ZW( N . M+ 1 ) . XW( N+ 1 , M+ 1  )  . 

1   YW(N+1 ,M+1 ) ,ZW(N+1 ,M+1) ,UTX(N.2) ,UTY(N.2) ,UTZ(N,2) .SX.SY.SZ.O) 
74  CONTINUE 

DO  76  N=1 ,NW 

■**  COMPUTE  INOUCED  VELOCITY  FROM  ONE  FORCESHOE  AT  THIS  RADIUS 

BUG=UTX(N.2)-UTX(N. 1)+U8X(N) 
CAT»UTY(N.2)-UTY(N. 1 )+UBY(N) 
D0G=UTZ(N.2)-UTZ(N, 1)+UBZ(N) 


UFP01540 
UFP01550 
UFP01560 
UFPO1570 
UFP01580 
UFP01590 
UFP01600 
UFP01610 
UFP01620 
UFP01630 
UFP01640 
UFP01650 
UFP01660 
UFP01670 
UFP01680 
UFP01690 
-UFP01700 
UFP01710 
UFP01720 
UFP01730 
UFP01740 
UFP01750 
UFP01760 
UFP01770 
UFP01780 
UFP01790 
UFP018CO 
UFP018 10 
UFP01820 
UFP01830 
UFP01840 
UFP01850 
UFP01860 
UFP01870 
UFP01880 
UFP01890 
UFP0190O 
UFP01910 
UFP01920 
UFP01930 
UFP01940 
UFP01950 
UFP01960 
UFP01970 
UFP01980 
UFP01990 
UFP020OO 
UFP02010 
UFP02020 
UFP02030 
UFP02040 
UFP02050 
UFP02060 
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IF(N.NE.NW)  GO  TO  77  UFP02070 

IF(M.NE.MID)  GO  TO  78  UFP02080 

C  UFP02090 

C  ***  ALL  TRANSITION  WAKE  VORTICES  DISAPPEAR  EXCEPT  MID,  WHICH  TURNS     UFP021OO 

C      INTO  ULTIMATE  WAKE  VORTICES  (TIP  AND  HUB)  BY  ASSUMING  MEAN  LOAD  I NGUFP02 1  10 

C  UFP02120 

BUG=BUG+WAKX+HBWX  UFP02  130 

CAT=CAT+WAKY+HBWY  UFP02  140 

D0G=D0G+WAK2+HBWZ  UFP02  150 

GO  TO  78  UFP02  160 

77  BUG=BUG-UBX(N+1 )  UFP02170 
CAT=CAT-UBY(N+1)  UFP02180 
D0G=00G-UBZ(N+1 )  UFP02190 

78  VOFWX(N,M)=BUG  ..  UFP0220O 
VOFWY(N.M)=CAT  UFP02210 
VOFWZ(N.M)=OOG  UFP02220 

76  CONTINUE  UFP02230 

92  DO  93  LN=1.NVFW  UFP02240 

KKM=K-LN  UFP02250 

IF  (KKM.LE.O)  KKM=KKM+NSR  UFP02260 

VIX=VIX+VOFWX(LN,M)-(GMEAN(M)  UFP02270 

1                   +OKFW(LN)*(GAM(M,KKM)-GMEAN(M) ) )  UFP02280 

VIY=VIY+VOFWY(LN.M)*(GMEAN(M)  UFP02290 

1                    +DKFW(LN)«(GAM(M.KKM)-GMEAN(M) ))  UFP02300 

93    VIZ=VIZ+VOFWZ(LN.M)»(GMEAN(M)  UFP02310 

1                    +DKFW(LN)*(GAM(M,KKM)-GMEAN(M) ) )  UFP02320 

71  CONTINUE  UFP02330 

'  10   RETURN  UFP02340 

END  UFP02350 


<<>' 
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SUBROUTINE 

RDFILE 

•    SUBROUTINE  RDFILE  READS  DISK  FILE14  OATA  CREATED  * 

•    BY  PUF- 

2  AND  XAVES  DATA  NEEDED  BY  FPUV  IN         * 

*    COMMONS 

GEOCOM  AND  SINCOM.   NOTE  THAT  NSR  IN      • 

*    PUF-2  MUST  BE  DIVISIBLE  BY  THE  NUMBER  OF  BLADES,  * 

*    NBLADE. 

*    NN 

=  >                                          * 

*    MM 

=  >                                          * 

•    NW 

=>  NUMBER  OF  VORTICES  IN  KEY  WAKE        * 

*    NBLADE 

=>  NUMBER  OF  PROPELLER  BLADES           * 

*    NVC 

=>  NO  OF  CHORD  VORTICES.  OTHER  BLADES   * 

•    NVS 

=>  NO  OF  SPAN  VORTICES.  OTHER  BLADES    * 

*    NVW 

=>  NO  OF  WAKE  VORTICES.  OTHER  WAKE      * 

*    NSR 

=>  NUMBER  OF  TIME  STEPS/REVOLUTION      * 

*    NWK 

=>  NUMBER  OF  RADII  IN  WAKE  HARMONICS    • 

•**  COMMON  TO  CONTAIN  GEOMETRIC  PARAMETERS 

COMMON  /  GEOCOM  /  MM . MN , NN , NW . NTORL . X ( 12 . 10) . Y ( 1 2 . 10) . Z ( 1 2 . 10 ) . 

5  XTIP(  10.  1 1 ).YTIP(  10.  11 ) .ZTIP(  10.  1 1 ) ,XH(  101 .7) ,YH(  101 .7). 

6  ZH( 101 .7) ,XW(21 . 10).YW(21 . 10) ,ZW(21 . 10) ,RH 

*»•  COMMON  TO  CONTAIN  SINGULARITY  STRENGTHS 

COMMON  /  SINCOM  /  GT ( 10 . 10 . 60 ) . G8 ( 10 . 9 . 60 ) . GTW ( 10 ) , SB ( 10 . 9 ) . GTV 
&   .GAM (9. 60) ,GMEAN(9) .DKFW(20) . H8WX . HBWY , HBWZ 

***  COMMON  TO  CONTAIN  FIELD  POINT  DATA.  NBLADE.  IDENT.  NSR.  AND  VSR 

COMMON  /  FPCOM  /  NBLADE . NSR . XFP . RFP . TZFP , YBAR , ZBAR . 
&    DBLAOE.DTFP. IDENT( 18) ,TFP(60) .VSR. IPLOT.DIAM 


DIMENSION  RZ( 10) ,R(9) .CH0CP(9) ,CHORZ( 10) .PHI B( 90) 

DIMENSION  X0(6.4.6).Y0(6,4,6).Z0(6.4.6) 

DIMENSION  RWK( 1 1 ) ,AWA( 16.11) ,BWA(  16.11) . AWR(  16,11)  . BWR (  16.  11) 
&   AWT( 16. 1 1 ) .BWT( 16. 1 1 ) 

DIMENSION  XW0(5.4.6) ,YW0(5.4,6) ,ZW0(5,4.6) 

DIMENSION  SB0(5.3) ,G8C(5.3.60) . GAMWC( 4 , 3 . 60) .GTC(5.4.60) . 
&   GAMC(3,60) ,GMEAN0(3) .DK0W(4) 

REA0(14)  AJ 

IF  (AJ.GE.2.0)  STOP 

RE AD ( 1 4 )  I STDY . NBLADE . MM . NN , NW . NVC , NVS . NVW , NSR , NWK 

•*»  SET  CONSTANT  FOR  GRID 


UFP02360 
UFP02370 
UFP02380 
UFP02390 
UFP0240O 
UFP02410 
UFP02420 
UFP02430 
UFP02440 
UFP02450 
UFP02460 
UFPO2470 
UFP02480 
UFP02490 
UFP025OO 
UFP02510 
UFP02520 
UFP02530 
UFP02540 
UFP02550 
UFP02560 
UFP02570 
UFP02580 
UFP02590 
UFP0260O 
UFP02610 
UFP02620 
UFP02630 
UFP02640 
UFP02650 
UFP02660 
UFP02670 
UFP02680 
UFP02690 
UFP027CO 
UFP02710 
UFP02720 
UFP02730 
UFP02740 
UFP02750 
UFP02760 
UFP02770 
UFP02780 
UFP02790 
UFP028CO 
UFP028  10 
UFP02820 
UFP02830 
UFP02840 
UFP02850 
UFP02860 
UFP02870 
UFP02830 
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NM=NN- 1 

MN=MM+1 

NP=NN+1 

NZ=NW+1 

NVCP=NVC+1 

NVSP=NVS+1 

NBL0=N8LA0E-1 

***  SET  CONSTANTS  FOR  UNSTEADY  PROBLEM 

IF  (ISTOY.EO. 1)  NSR=1 

NVVP=NVW+1 

NTB=NN-MM 

***  READ  KEY  BLADE  GEOMETRY 


M=1  .MN), 
M=1  .MN) , 

MN)  . 


REAO(14)  ( (X(N,M) ,N=1 ,NP) 

2  ( (Y(N.M) ,N=1 ,NP) 

3  ( (Z(N.M) ,N»1 ,NP) ,M= 1 

4  •  (RZ(M) ,M=1 .MN), 

5  (R(M) ,M»1 ,MM) , 

6  (CHOCP(M) , M= 1 .MM) , 

7  (CHORZ(M) ,M= 1 ,MN) , 

8  (PHIB(K) ,K= 1 .NTB) 

■•  READ  BLADE  GEOMETRY  OF  OTHER  BLADES 


IF  (NBLADE.LE. 1 )  GOTO 
READ(14)  ( ( (XO(N.M.K) , 

2  (((YO(N.M.K) 

3  ( ((ZO(N.M.K), 
11   CONTINUE 


1  1 

N=1  .NVCP) ,M=1 .NVSP) 
N=1  .NVCP) ,M=1 .NVSP) 
N=1  .NVCP) .M=1  .NVSP) 


K=1  .N8L0) 
K=  1  .N8L0) 
K=1  .N8L0) 


»**  READ  WAKE  RADII  AND  WAKE  HARMONICS 


READ(14)  (RWK(M) ,M=1 . NWK ) 


READ(14)  ( (AWA(I .M) .1= 

2  ( (BWA(I.M). 1= 

3  ( (AWR(I .M) . 1= 

4  ( (BWR( I .M) . 1= 

5  ((AWT(I ,M),I= 

6  ((BWT(I ,M),I= 


.  16).M=1 .NWK) 
.  16) .M=1 .NWK) 
.  16) ,M=1 ,NWK) 
.  16)  .M=1 .NWK) 
.  16)  ,M=1 ,NWK) 
.  16)  .M=1 ,NWK) 


*»*  READ  VARIOUS  PROPELLER/WAKE  CONSTANTS 

READ( 14)  WAKE.RPM.DIAM.RH.CDRAG.VSR.UR. ( IDENT(N) . N=1  .  18) 


•»*  REAO  KEY  WAKE  GEOMETRY 


REAO(  14) 


((XW(N.M) 
((YW(N.M) 
((ZW(N.M) 


N=1  .NZ)  ,M=»1 
N=1 ,NZ) ,M=1 
N»1,N2).M»1 


MN)  . 
MN)  . 
MN)  . 


( (XH(N.K) ,N»1 . 101 ) ,K=1 .N8LADE) . 


UFP02890 
UFP02900 
UFP02910 
UFP02920 
UFP02930 
UFP02940 
UFP02950 
UFP02960 
UFP02970 
UFP02980 
UFP02990 
UFP030O0 
UFP03010 
UFP03020 
UFP03030 
UFP03040 
UFP03050 
UFP03060 
UFP03070 
UFP03080 
UFP03090 
UFP031CO 
UFP031 10 
UFP03120 
UFP03130 
UFP03140 
UFP03150 
UFP03  160 
UFP03170 
UFP03180 
UFP03190 
UFP0320O 
UFP03210 
UFP03220 
UFP03230 
UFP03240 
UFP03250 
UFP03260 
UFP03270 
UFP03280 
UFP03290 
UFP0330O 
UFP03310' 
UFP03320 
UFP03330 
UFP03340 
UFP03350 
UFP03360 
UFP03370 
UFP03380 
UFP03390 
UFP0340O 
UFP034  10 
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5  ((YH(N.K).N=1. 101 ),K=1 .N8LA0E) .  UFP03420 

6  ( (ZH(N.K) ,N=1 , 101 ) ,K=1 .NBLADE) .  UFP03430 

7  ((XTIP(N.L).L=1 ,NP) .N=1 .NN)  .  UFP03440 

8  ((YTIP(N.L) ,L=1 ,NP) ,N=1 ,NN) ,  UFP03450 

9  ( (ZTIP(N.L) ,L=1 ,NP) ,N=1 .NN)  ,  UFP03460 
C  UFP03470 
C  ••»  READ  KEY  8LADE  SINGULARITY  STRENGTHS  UFP03480 
C  UFP03490 

•  ( (SB(N.M) .N=1 ,NN) ,M=1 , MM ) .  UFP03500 

1  ( ( (GB(N.M.KK) ,N=1 .NN) ,M=1 .MM) ,KK=1 .NSR) .  UFP03510 

2  ( ( (GT(N.M.KK) .N=1 ,NN) ,M=1 ,MN) .KK=1 .NSR).  UFP03520 

3  ((GAM(M.KK) ,M=1 .MM) .KK=1 ,NSR)  UFP03530 
C  UFP03540 
C  »**  OTHER  WAKE  GEOMETRY  ..  UFP03550 
C  UFP03560 

IF  (NBLADE. LE.1)  GOTO  12  UFP03570 

READ(14)  ( ( (XWO(N.M.K) ,N*1 ,NVWP) .  M=  1 .NVSP) ,K= 1 .NBLO) .  UFP03580 

2  ( ((YWO(N.M.K) ,N=1  .NVWP) ,M=1 .NVSP) ,K=1  .NBLO)  ,  UFP03590 

3  ( ((ZW0(N.M.K).N=1 .NVWP) ,M=1 .NVSP) .K=1 .NBLO) .  UFP0360O 
C  .  UFP036  10 
C  »**  OTHER  BLADE  SINGULARITIES  UFP03620 
C  UFP03630 

4  (  (SBO(N.M)  ,N=1  .NVC)  ,M=<1  .NVS)  ,  UFP03640 

5  ( ((GBC(N.M.KK) ,N=1 .NVC) .M=1 ,NVS) ,KK=1 .NSR) .  UFP03650 

6  ( ( (GAMWC(N.M.KK) .N=1 .NVW) ,M=1 ,NVS) .KK=1 .NSR) .  UFP03660 

7  (( (GTC(N.M.KK) ,N=1 .NVC) ,M=1 .NVSP) ,KK=1 , NSR) .  UFP03670 

8  ( (GAMC(M.KK) ,M=1 .NVS) ,KK=1 .NSR) ,  UFP03680 

9  (GMEAN(M) ,M=1 .MM) ,  UFP03690 

*  (GMEANO(M) ,M=1 ,NVS) .  UFP03700 

1  (DKFW(N) ,N=1 ,NW) .  UFP03710 

2  (DKOW(N) .N=1 .NVW)  .  UFP03720 
12  CONTINUE  UFP03730 
99   GOTO  100  UFP03740 

C  UFP03750 

C  &&&  DIAGNOSTICS  &&&&&&&&&&&&&&&&&&&&&&&    ,  UFP03760 

C  UFP03770 

901  F0RMAT(16I5)  UFP03780 

902  F0RMAT(8F10.5)  UFP03790 

903  F0RMAT(18A4)  UFP038OO 
WRITE (1,902)  AJ  UFP03810 
WRITE(1,901)  ISTDY. NBLADE .MM. NN , NW . NVC. NVS. NVW. NSR, NWK  UFP03820 
WRITE(1.902)  ((X(N.M) ,N=1 . NP ) ,M=1 ,MN).  UFP03830 

2  ( (Y(N.M) .N=1 .NP) ,M=1 ,MN) .  UFPO3840 

3  ((Z(N.M) ,N=1 ,NP) ,M=1 .MN).  UFP03850 

4  (RZ(M) ,M»1 .MN) ,  UFP03860 

5  (R(M) ,M=1 .MM) .  UFP03870 

6  (CHOCP(M) ,M=1 .MM) ,  UFP03880 

7  (CHORZ(M) ,M=1 .MN) .  UFP03890 

8  (PHIB(K) ,K»1 ,NTB)  UFP03900 
IF  (N8LADE.LE. 1)  GOTO  21  UFP03910 
WRITE ( 1 .902)  (( (X0(N,M.K).N=1 .NVCP) . M« 1 .NVSP) ,K=1 .NBLO) .  UFP03920 

2  (((YO(N.M.K) .N=1 .NVCP) ,M=1 .NVSP) ,K=1 .N8L0) .  UFP03930 

3  (((ZO(N.M.K) ,N=1 .NVCP) .M=1 .NVSP) .K«1 .NBLO)  UFP03940 
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21 


22 
100 


CONTINUE 

WRITE(1.902)  (RWK(M) ,M=1 , NWK ) 

WRITE(1.902)  ((AWA(I.M) 


WRITE(  1 
WRITE(  1 
WRITE(  1 


( (BWA( I ,M) . 1=1 
■  ((AWR(I ,M) . 1=1 
((BWR(I.M).I=1 
( (AWT(I.M) ,1=1 
( (BWT(I.M). 1=1 

902)  WAKE.RPM.DIAM 

903)  (IDENT(N) ,N=1 
902)  ( (XW(N.M) ,N=1 


1=1 
16) 
16) 
16) 
16) 
16) 


. 16) ,M=1 .NWK) , 

,M=1 .NWK) . 

.M=1 .NWK) , 

,M=1 .NWK) . 

,M= 1 .NWK ) . 

,M=1 .NWK) 

RH.CDRAG.VSR.UR 

18) 

NZ)  ,M=1 .MN)  . 


2  ( (YW(N.M) ,N=1 ,NZ) ,M=1 .MN) . 

3  ( (ZW(N,M).N=1 .NZ) ,M=1 . MN ) . 

4  (  (XH(N.K),N*1  .  10.1  )  ,K=1  .NBLAOE)  . 

5  ((YH(N.K) ,N=1, 101).K=1 .NBLADE) . 
.6  ((ZH(N.K).N=1 . 101 ) .K=1 .NBLADE)  . 

7  ((XTIP(N.L) .  L=1 ,NP).N=1,NN), 

3  ((YTIP(N.L) ,L=1 ,NP) ,N=1 .NN). 

9  ((ZTIP(N.L) ,L=1 ,NP) ,N=1 ,NN) 

WRITE ( 1 ,902)  ( (SB(N.M) .N=1 .NN) ,M=1 .MM) . 

2  ( ( (GB(N.M.KK) ,N=1 .NN) ,M=1 .MM) .KK=1 .NSR)  . 

3  ( ( (GT(N.M.KK) ,N=1 , NN ) ,M=1 ,  MN  ) ,KK=1 .NSR) . 

4  ( (GAM(M.KK) ,M=1 .MM) ,KK=1 ,NSR) 
IF  (NBLADE. LE.  1  )  GOTO  22 

WRITE (  1 .902)  ( ( (XWO(N.M.K) ,N= 1  , NVWP ) ,M=1  ,NVSP)  , K=1  ,NBLO)  . 

2  .    ( ( (YWO(N.M.K) ,N=1  ,NVWP) ,M=1 ,NVSP) ,K=1 ,NBLO)  , 

3  (( (ZWO(N.M.K) ,N=1 .NVWP) ,M=1 .NVSP) ,K=1 .NBLO) 
WRITE (  1 .902)  ( (SBO(N.M) ,N=1 .NVC) ,M=1 ,NVS)  . 


CONTINUE 

RETURN 

ENO 


( ( (GBC(N.M.KK) ,N= 1 .NVC) ,M=1 . NVS ) ,KK=1 ,NSR) . 
( ( (GAMWC(N.M.KK) ,N=1 . NVW ) ,M=1 .NVS) .KK=1 ,NSR) . 
( ( (GTC(N.M.KK) ,N=1 , NVC ) ,M=1 ,NVSP) ,KK=1 .NSR) . 
( (GAMC(M.KK) ,M=1 .NVS) ,KK=1 ,NSR) , 
(GMEAN(M) ,M=1 .MM) . 
( GME ANO ( M ) , M= 1 . NVS )  , 
(DKFW(N).N=1 ,NW) . 
(DK0W(N).N=1 ,NVW) 


/ 


UFP03950 
UFP03960 
UFP03970 
UFP03980 
UFP03990 
UFP040O0 
UFP04010 
UFP04020 
UFP04030 
UFP04040 
UFP04050 
UFP04060 
UFP04070 
UFP04080 
UFP04090 
UFP04 100 
UFP04 1 10 
UFP04 120 
UFP04 130 
UFP04140 
UFP04 150 
UFP04 160 
UFP04170 
UFP04 180 
UFP04 190 
UFP0420O 
UFP04210 
UFP04220 
UFP04230 
UFP04240 
UFP04250 
UFP04260 
UFP04270 
UFP04280 
UFP04290 
UFP0430O 
UFP04310 
UFP04320 
UFP04330 
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SUBROUTINE  RDTERM( ISTOP) 


C 
C 
C 

c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


901 
10 


20 
902 


SUBROUTINE  RDTERM  READS  IN  BLADE  AND  GEOMETRY 
DATA  INTERACTIVELY  AND  COMPUTES  NECESSARY  FIELD 
POINT  GEOMETRY. 

NBLADE       =>  NUMBER  OF  BLAOES 

NTORL   ■  0  =>  CONSIDER  THICKNESS  ANO  LOADING 

*  1  =>  CONSIDER  THICKNESS  ONLY 

«  2  =>  CONSIDER  LOADING  ONLY 
XFP         =>  AXIAL  FIELD  POINT  COORDINATE. 

POSITIVE  OOWNSTREAM 
RFP         =>  RAOIAL  flELD  POINT  COORDINATE 
TZFP        =>  INITAL  FIELD  POINT  ANGLE  IN  DEG. 
Y8AR.ZBAR   =>  Y  AND  Z  FIELD  POINT  COORDINATES 
TH  =>  ANGLE  FROM  TIME  STEP  INPUT 

YB.ZB        =>  Y.Z  FIELD  PT  COORD.  BEFORE  TIME  ST 


'*  COMMON  TO  CONTAIN  GEOMETRIC  PARAMETERS 

COMMON  /  GEOCOM  /  MM . MN , NN . NW . NTORL . X (  1  2 .  10 ) . Y (  1 2 .  1 0) , 2 (  1 2 .  10 ) . 
&    XTIP(  10.  1  1  )  . YTIP(  10.  1  1 ) .ZTIP(  10.  1 1 ) .XH(  101 .7) . YH(  101 .7). 
&    ZH( 101 .7) ,XW(21 . 10) .YW(2  1 . 10) .ZW(21 . 10) ,RH 

-  COMMON  TO  CONTAIN  FIELD  POINT  DATA,  NBLAOE .  IDENT.  NSR .  ANO  VSR 

COMMON  /  FPCOM  /  NBLADE . NSR . XFP . RFP . TZFP , YBAR . ZBAR , 
&    DBLADE.DTFP. IDENT (  18) ,TFP(60) .VSR. I  PLOT. D I  AM 

FORMAT(/'  ENTER  XFP,  RFP,  TZFP.'/ 
&    '  FOR  HELP.  ENTER  -11.  FOR  XFP.') 
WRITE(6.901) 
READ(5,»)  XFP, RFP, TZFP 
IF  (XFP.LT. -10. )  GOTO  90 
IF  (TZFP. LT. 0.0)  GOTO  70 
IF  (RFP.GT.O. )  GOTO  60 
WRITE(6,902) 

FORMAT(/'  ENTER  NBW , N , M , KK ' / '  FOR  HELP.  ENTER  ZERO  FOR  NBW.') 
REA0(5.«)  NBW.N.M.KK 

TH=6. 283 185* FLOAT (KK-1 )/FLOAT(NSR) 
IF  (NBW.EO.O)  GOTO  80 

■»  FIELD  POINT  CENTERED  ON  BLADE  GRID. 

IF(NBW.EO.I)  CALL  FPGRID( X , Y . Z . N , M . KK . 12 ) 
■•  FIELD  POINT  CENTERED  ON  TRANSITION  WAKE  GRID. 

IF(N8W.E0.2)  CALL  FPGRID(XW , YW, ZV , N. M, KK , 2 1 ) 


UFP04340 
UFP04350 
UFP04360 
UFP04370 
UFP04380 
UFP04390 
UFP0440O 
UFP04410 
UFP04420 
UFP04430 
UFP04440 
UFP04450 
UFP04460 
UFP04470 
UFP04480 
UFP04490 
UFP045OO 
UFP04510 
UFP04520 
UFP04530 
UFP04540 
UFP04550 
UFP04560 
UFP04570 
U.rP04580 
UFP04590 
UFP046O0 
UFP04610 
UFP04620 
UFP04630 
UFP04640 
UFP04650 
UFP04660 
UFP04670 
UFP04680 
UFP04690 
UFP0470O 
UFP04710 
UFP04720 
UFP04730 
UFP04740 
UFP04750 
UFP04760 
UFP04770 
UFP04780 
UFP04790 
UFP048OO 
UFP04810 
UFP04820 
UFP04830 
UFP04840 
UFP04850 
UFP04860 
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IF(N8W.NE. 1 . AND.NBW.NE.2)  GOTO  10 
60   DTFP=36O.O/FL0AT(NSR) 

RETURN 
70   IST0P=1 

RETURN 
80   WRITE(6.903) 

903  FORMAT(/'  NEW  =  0  =>  HELP'/5X. 
&   '=  2  =>  TRANSITION  WAKE  GRID 
&   '  M',7X.'=>  SPANWISE  COORDINATE'/' 
&   '=>  TIME  STEPS  ROTATED  +  1'/) 

GOTO  20 
90   WRITE(6,904) 

904  FORMAT(/'  XFP     <-10=>  HELP'/8X. 
AXIAL  FIELD  POINT  COORDINATE.  POSITIVE  DOWNSTREAM'/ 

<  0  =>  SELECT  COORDINATE  USING  BLADE  OR  WAKE  GRID'/ 
=>  RADIAL  FIELD  POINT  COORDINATE'/ 

<  0  =>  ST0P'/8X.'>  0  =>  ANGULAR  FILE  POINT  COORDINATE 


■  1  =>  BLADE  WAKE  GRID'/5X, 
/'  N'.7X,'=>  CHOROWISE  COORDINATE 
KK' ,6X. 


'>- 10=> 
'  RFP 
8X. '>  0 
'  TZFP 
/) 

GOTO  10 

END 


UFP04870 
UFP04880 
UFP04890 
UFP049OO 
UFP04910 
UFP04920 
UFP04930 
/UFP04940 
UFP04950 
UFP04960 
UFP04970 
UFPO4930 
UFP04990 
UFP05000 
UFP05010 
UFP0502O 
UFP05030 
UFP05040 
UFP05050 
UFP05060 
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SUBROUTINE  FPGRID(XX.YY.ZZ.N.M,KK,NDIM) 


SUBROUTINE  GRID  CALCULATES  X.R.  AND  THETA  FIELD  • 

POINT  GEOMETRY  COMPONENTS  GIVEN  GRID  VECTORS  * 

X,Y,  AND  Z  FOR  BLADE  OR  WAKE,  GRID  COORDINATES  * 

N  AND  M,  AND  TIME  STEP  KK .  • 

■ 
XX.YY.ZZ   =>  GRID  GEOMETRY  VECTORS 

N.M        =>  GRID  INDICES  SELECTED  FOR  VELOCITY  * 

KK          =>  TIME  STEP.  NUMBER  OF  ANGLE  STEPS  +  1* 

XFP, Y8AR. ZBAR.TZFP.RFP  * 

=>  X.Y.Z. THETA. RADIAL  FIELD  PT  COORD  * 

DBLAOE     =>  ANGLE  BETWEEN  BLADES  * 

TH         =>  ANGLE  OF  TIME  STEP  * 

NDIM       =>  FIRST  DIMENSION  OF  GRIDS  PASSED  • 


C 
C 
C 

c 

C 
C 
C 
C 

c 
c 
c 
c 

C   '   ■ 

c 

C 

c 

C      .ft..*,..,......,,...,,...,......,...............*...., 

C 

C  ***  COMMON  TO  CONTAIN  FIELD  POINT  DATA,  NBLADE.  IDENT.  NSR .  AND  VSR 
C 

COMMON  /  FPCOM  /  NBLADE , NSR . XFP , RFP . TZFP , YBAR , ZBAR , 
&    D8LADE,DTFP.IDENT( 18) , TFP ( 60) . VSR , IPLOT , DI AM  . 
C 

DIMENSION  XX (NDIM, 10) , YY(NDIM, 10) ,ZZ(NDIM, 10) 
C 

C  *»«  COMPUTE  X.Y.Z  COORDINATES  IN  ZERO  ANGLE  POSITION 
C 

TH=6.283  185«FL0AT(KK-1 ) /FLOAT (NSR) 

THDEG=TH*360./6.283  185 

DBLADE =360.0/ FLOAT (NBLADE) 

XFP=0.25*(XX(N.M)  +  XX(N+1.M)  +  XX(N.M+1)  +  XX(N+1.M+1)) 

YB  =0.25*(YY(N,M)  +  YY(N+1.M)  +  YY(N.M+1)  +  YY(N+1,M+1)) 

ZB  =0.25*(ZZ(N.M)  +  ZZ(N+1,M)  +  ZZ(N,M+1)  +  ZZ(N+1.M+1)) 
C 

C  ***  CONVERT  TO  CYLINDRICAL  COORDINATES 
C 

RFP=S0RT(YB*Y8+ZB*ZB) 

TZFP=ATAN2(ZB.YB)/1 .7453293E-02 

TZFP=TZFP+THOEG 
C 

C  ***  BRANCH  IF  GENERATING  PIEWAKE  DATA 
C 

IF  (IPLOT. EO. 2)  GOTO  20O 

TZFP=180.0/FL0AT(N8LADE)  -  ATAN2 ( ZB . YB )/ 1 . 7453293E -02 
C 

C  **•  KEEP  TZFP  BETWEEN  0  AND  ONE  BLADE  ANGLE 
C 

DO  100  K=1 .NBLADE 

IF(TZFP.LT.O.O)  TZFP=TZFP+DBLADE 

IF(TZFP.GT. DBLADE)  TZFP-TZFP -DBLADE 
100   CONTINUE 
C 


UFP05070 
UFP05080 
UFP05090 
UFP0510O 
UFP051 10 
UFP05120 
UFP05130 
UFP05140 
UFP05150 
UFP05160 
■ UFP05170 
UFP05180 
UFP05190 
UFP05200 
UFP05210 
UFP05220 
UFP05230 
UFP05240 
UFP05250 
UFP05260 
UFPO5270 
UFP05280 
UFP05290 
•  UFP05300 
UFP05310 
UFP05320 
UFP05330 
UFP05340 
UFP05350 
UFP05360 
UFPO5370 
UFP05380 
UFP05390 
UFP05400 
UFP05410 
UFP05420 
UFP05430 
UFP05440 
UFP05450 
UFP05460 
UFP05470 
UFP05480 
UFP05490 
UFP05500 
UFP05510 
UFP05520 
UFP05530 
UFP05540 
UFP05550 
UFP05560 
UFP05570 
UFP05580 
UFP05590 
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C  *»*  ROTATE  BLADE  FOR  PIEWAKE  DATA  UFP05600 

C  UFPC56  10 

200   YBAR=Y8*C0S(TH)  -  ZB*SIN(TH)  UFP05620 

ZBAR=YB»SIN(TH)  +  ZB-COS(TH)  UFP05630 

30O   RETURN  UFP05640 

END  UFP05650 


/ 
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c 
c 

c 
c 
c 
c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


SUBROUTINE  FPSTEP(K) 


.•••-............ 


FPSTEP  IS  A  SUBROUTINE  TO  CALCULATE  THE  FIELD 
POINT  VELOCITIES  INDUCED  BY  ONE  BLADE  AT  "TIME 
STEP*.  K.   THE  "DIRTY  WORK"  IS  DONE  IN 
SUBROUTINE  FPVEL. 

TFP       =>  ANGULAR  FIELD  POINT  COORDINATE.  DEC 
TRAD      =>  ANGULAR  FIELD  POINT  COORDINATE.  RAD. 
YFP.2FP   =>  Y  AND  Z  FIELD  POINT  COORDINATES 
VIX.VIY.VIZ.VIR.VIT  =>  INOUCED  VELOCITY  COMPONENT 


»*»  COMMON  TO  CONTAIN  GEOMETRIC  PARAMETERS 

COMMON  /  GEOCOM  /  MM . MN. NN . NW , NTORL , X (  1 2 .  10) . Y ( 1 2 .  10 ) . Z (  1  2 .  10  )  . 
&    XTIP(10.11).YTIP(10.11).ZTIP(10.11).XH(101,7) ,YH( 101 .7). 
&    ZH( 101 ,7).XW(21 . 10) .YW(21 , 10) ,ZW(21 . 10).RH 

*•»  COMMON  TO  CONTAIN  SINGULARITY  STRENGTHS 

COMMON  /  SINCOM  /  GT ( 10. 10 , 60) . G8 ( 10 . 9 . 60 ) , GTW( 10 ) . SB ( 10 . 9 ) . GTV 
&   .GAM (9. 60) ,GMEAN(9) .  DKFW(20) . HBWX . HBWY . HBWZ 

***  COMMON  TO  CONTAIN  FIELD  POINT  DATA.  N8LADE.  IDENT.  NSR .  ANO  VSR 

COMMON  /  FPCOM  /  NBLADE.NSR.XFP.RFP.TZFP. YBAR, ZBAR. 
&    DBLADE.DTFP, IDENT ( 18) . TFP( 60) . VSR . I  PLOT. DI AM 

»**  COMMON  TO  CONTAIN  VELOCITY  DATA 

COMMON  /  VELCOM  /  VIX( 60 ) . VIT ( 60 ) , VI R ( 60) . U( 60. 3 ) . NSAMP 

»•*  CALCULATE  CARTESIAN  FIELD  POINT  COORDINATES  TO  FEED  FPVEL 

TFP(K)=FL0AT(K-1)-DTFP-TZFP 

TRAD=TFP(K)*1.7453293E-02 

COSFP=COS(TRAD) 

SINFP*SIN(TRAD) 

YFP=RFP*COSFP 

ZFP=RFP*SINFP 

•*•  FPVEL  WILL  CALCULATE  VIX,  VIY.  ANO  VIZ  FOR  FIELD  POINT 

CALL  FPVEL(K,XFP,YFP.ZFP.VIX(K) .VIY. VIZ) 

VIR(K)=VIY«COSFP  +  VIZ'SINFP 

VIT(K)=VIZ*COSFP  -  VIY«SINFP 

RETURN 

END 


UFP05660 
UFP05670 
UFP05680 
UFP05690 
UFP05700 
UFP057  10 
UFP05720 
UFP05730 
UFP05740 
UFP05750 
UFP05760 
UFP05770 
UFP05780 
UFP05790 
UFP0580O 
UFP05810 
UFP05820 
UFP05830 
UFP05840 
UFP05850 
UFP05860 
UFP05870 
UFP05880 
UFP05890 
UFP0590O 
UFP05910 
UFP05920 
UFP05930 
UFP05940 
UFP05950 
UFP05960 
UFP05970 
UFP05980 
UFP05990 
UFP060OO 
UFP06010 
UFP06020 
UFP06030 
UFP06040 
UFP06050 
UFP06060 
UFP06070 
UFP06080 
UFP06090 
UFP06100 
UFP06  1 10 
UFP06  120 
UFP06130 
UFP06  140 
UFP06  150 
UFP06160 
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SUBROUTINE  SUMVEL  UFP06170 

C      ,,.....,«.«...,.....«.............»....,,..,....«...««  UFP06  180 

C      •  *              UFP06190 

C      •    SUBROUTINE  SUMVEL  SUMS  VELOCITIES  FROM  EACH  *              UFP0620O 

C      *    BLADE.  CALCULATED  IN  FPSTEP.  TO  PRODUCE  THE  •              UFP062  10 

C      *    TOTAL  PROPELLER-INOUCED  FIELD  POINT  VELOCITY.  *              UFP06220 

C  *              UFP06230 

C      *    VIX.VIR.VIT  =>  VELOCITY  COMPONENTS  INDUCED  BY  »              UFPO6240 

C      •                     ONE  BLADE  AT  "TIME  STEP*  K  *              UFP06250 

C      *    U(N.d)       ■>  VELOCITY  INDUCED  BY  ALL  BLAOES  •              UFP06260 

C      *                     AT  'TIME  STEP"  N  •              UFP06270 

C    '  *            J=0  =>  AXIAL  COMPONENT  «              UFP06280 

C      *             J-1  =>  RADIAL  COMPONENT  *              UFP06290 

C      *             J=2  =>  TANGENTIAL  COMPONENT  •              UFP063CO 

C      *    NSAMP         =>  NUMBER  OF  'TIME  STEPS"  IN  ONE  *              UFP063  10 

C      »                    BLADE  ANGLE  (=  NUM8ER  OF  *              UFP06320 

C      »                    U  VELOCITIES  COMPUTED.)  *              UFP06330 

C      •  *              UFP06340 

C      ......................................................  UFP06  35O 

C  UFP06360 

C  »-*  COMMON  TO  CONTAIN  FIELD  POINT  DATA.  NBLADE .  IDENT.  NSR.  AND  VSR    UFP06370 

C  UFPO6330 

COMMON  /  FPCOM  /  NBLADE . NSR ,XFP . RFP , TZFP , YBAR  .  ZBAR  ,  UFP06390 

&    DBLADE.DTFP. IDENT( 18),TFP(60) .VSR. IPLOT.DIAM  UFP0640O 

C  UFP064  10 

C  *»*  COMMON  TO  CONTAIN  VELOCITY  DATA                       .  UFP06420 

C  UFP06430 

COMMON  /  VELCOM  /  VIX ( 60 ) . VI T ( 60) . VI R( 60 ) . U( 60 . 3 ) . NS AMP  UFPO6440 

C  UFP06450 

C  **-  AT  ANY  GIVEN  REAL  TIME.  EACH  BLADE  CAN  BE  CONSIDERED  TO  HAVE        UFP06460 

C      TRAVELLED  THROUGH  A  CERTAIN  NUMBER  OF  "TIME  STEPS".  UFP06470 

C    '   BY  SUMMING  VIX.  VIR.  AND  VIT  AT  8LADE  RATE.  WE  CAN  INCLUDE          UFP06480 

C      THE  VELOCITIES  INDUCED  BY  EACH  BLADE  INTO  ONE  VELOCITY  UFP06490 

C      INDUCED  BY  THE  ENTIRE  PROPELLER.  UFP065OO 

C                  • •  UFP065  10 

NSAMP=NSR/NBLADE  UFP06520 

DO  20  N=1. NSAMP  UFP06530 

DO  10  0=1,3  UFP06540 

10   U(N,U)=0.0  UFP06550 

DO  20  K=1. NBLADE  UFP06560 

L=N+(K-1 )*NSAMP  UFP06570 

U(N. 1)=U(N.  1  )+VIX(L)  UFP06580 

U(N,2)=U(N,2)+VIR(L)  UFP06590 

20   U(N.3)=U(N,3)+VIT(L)  UFP066CO 

C  UFP066  10 

C  *«•  NOW  CONVERT  FROM  REFERENCE  VELOCITY  TO  SHIP  VELOCITY  UFP06620 

C  UFP06630 

DO  30  N=1. NSAMP  UFP06640 

DO  30  J=1 .3  UFP06650 

30   U(N. J)»U(N. J)/VSR  UFP06660 

RETURN  UFP06670 

ENO  UFP06680 
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SUBROUTINE  FOUCAL 


«    SUBROUTINE 

FOUCAL  TAKES  THE  PROPELLER  INOUCED 

*    VELOCITIES 

CALCULATED  IN  SUMVEL  AND  PERFORMS 

»    A  HARMONIC 

ANALYSIS  ON  THEM.   IT  CALCULATES        * 

»    FOURIER  COEFFICIENTS  FOR  BOTH  COSINE  (A)  AND      * 

*    SINE  (B)  AND  ALSO  COSINE  AMPLITUDE  (AMP)  AND      * 

*    PHASE  ANGLE  (PH) 

*    NH 

=  > 

NUMBER  OF  HARMONICS  CALCULATED   * 

*    A(N,J).B(N, 

J)  => 

FOURIER  COEFFICIENTS              * 
A  =>  COSINE  COEFFICIENTS 
B  =>  SINE  COEFFICIENTS           - 
N  =>  HARMONIC  NUMBER  UP  TO  15    * 

J  = 

1  =>  AXIAL  VELOCITY 

3 

2  =>  RADIAL  VELOCITY              • 

3 

3  =>  TANGENTIAL  VELOCITY 

*    AMP (N.J) 

=  > 

FOURIER  COEFFICIENT  AMPLITUDE    * 
N  AND  J  SUBSCRIPTS  AS  FOR  A,  B   * 

•    PH(N.J) 

=  > 

PHASE  ANGLE  ASSOCIATED  WITH  AMP  * 

*    NFB 

=  > 

NUMBER  OF  FIELD  POINTS  /  BLAOE   * 

»**  COMMON  TO  CONTAIN  FIELD  POINT  DATA.  NBLADE.  IDENT,  NSR,  AND  VSR 


C 

C  »- 

c 
c 

c  ** 

c 
c 

c  •* 

c 


10 


COMMON  /  FPCOM  /  NBLADE , NSR , XFP , RFP , T2FP . YBAR , 2BAR , 
&    DBLADE.DTFP. IDENT ( 18).TFP(60) .VSR. I  PLOT. D I  AM 

*  COMMON  TO  CONTAIN  VELOCITY  DATA 

COMMON  /  VELCOM  /  VIX ( 60 ) . VIT( 60) . VIR( 60) . U( 60 . 3 ) . NSAMP 

•  COMMON  TO  CONTAIN  HARMONICS  COEFFICIENTS 

COMMON  /  HARCOM  /  A(  1 5 . 3  )  . B (  1  5 . 3  )  . AMP (  1 5 . 3 ) . PH(  1 5 . 3 ) . NH 

-  FIRST  CALCULATE  A  ANO  B 

NFB=NSR/NBLAOE 

NHMAX=NFB/2+1 

NH=MINO(NHMAX. 15) 

STEP=6.283185/NFB 

ANFB=2.0/NFB 

BNFB=1 .O/NFB 

DO  11  K=1,NH 

DO  10  J=1 .3 

A(K.J)=0.0 

B(K.J)=0.0 

DO  9  N»1 .NFB 

T=((N-1)*(K-1))-STEP 

ST=SIN(T) 


UFP06690 
UFP067CO 
UFPO6710 
UFP06720 
UFPO6730 
UFP06740 
UFP06750 
UFP06760 
UFP06770 
UFP06780 
UFP06790 
UFP0680O 
UFP06810 
UFP06820 
UFP06830 
UFP06840 
UFP06850 
UFP06860 
UFP06870 
UFP06880 
UFP06S90 
UFP069CO 
UFP059 10 
UFPO6920 
UFP06930 
UFP06940 
UFP06950 
UFP06960 
UFP06970 
UFP06980 
UFP06990 
UFP07000 
UFP07010 
UFP07020 
UFP07030 
UFP07040 
UFP07050 
UFP07060 
UFP07O70 
UFP07080 
UFP07090 
UFP07  100 
UFP07  1 10 
UFP07120 
UFP07130 
UFP07140 
UFP07150 
UFP07160 
UFP07170 
UFP07180 
UFP07 190 
UFP0720O 
UFP07210 


-  80  - 


CT=COS(T)  UFP07220 

DO  9  J=1  .3  UFP07230 

A(K. J)=A(K,J)+U(N, J)*CT  UFP07240 

9     B(K. J)=B(K, J)+U(N. J)*ST  UFP07250 

C  UFP07260 

C  •»*  CONVERT  A  AND  B  INTO  AMP  AND  PH  UFP07270 

C  UFP07280 

DO  11  J=1,3  UFP07290 

IF(K.E0.1)  A(  1 . J)=A(  1  , J)»BNFB  UFP073OO 

IF(K.EO.NHMAX)  GO  TO  12  UFP07310 

IF(K.GT.1)  A(K. J)=A(K.d)*ANFB  UFP07320 

IF(K.GT.1)  B(K. J)=B(K, J)»ANFB  UFP07330 

GO  TO  13  UFP0734O 

12  A(K,d)=A(K,  J)*BNFB  ..  UFP07350 
B(K.J)=0-0  UFP07360 

13  AMP(K, J)=SQRT(A(K, J)«-2+B(K, J)**2)  UFP07370 
PH(K, J)=ATAN2(A(K,0) , B ( K , J ) ) *57  .  296  UFP07380 
IF(PH(K.J) .LT.O.O)  PH(K,J)=PH(K, J)+360.0  UFP07390 
PH(K.  J)=>PH(K,  J)/NBLADE  UFP074OO 
PH(1.J)«0.0  UFP07410 
AMP(  1  . J)»A(  1 . J)  UFP07420 

11    CONTINUE  UFP07430 

RETURN  UFP07440 

END  UFP07450 


y 
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SUBROUTINE  PI EDAT( ISTOP ) 


SUBROUTINE  PIEDAT  GENERATES  DATA  FROM  INPUT  TO 
AXIAL  VELOCITY  OUTPUT  TO  BE  USED  AS  INPUT  FOR 
UNSTEADY  PIEWAKE.   IT  USES  BLADE  GRID  GEOMETRY 
TO  CALCULATE  NSR  FIELD  POINT  AXIAL  T IME -  AVERAGED 
VELOCITIES  AT  FIVE  RADII  ALONG  THE  LEADING  EDGE. 
IT  THEN  CALCULATES  NSR  FIELD  POINT  AXIAL  TIME- 
AVERAGED  VELOCITIES  AT  1.1  TIMES  PROPELLER  RADIUS 

N  =>  CHORDWISE  GRID  COORDINATE 

M  =>  SPANWISE  GRID  COORDINATE 

K.KK        =>  ANGLE  INCREMENTS  +  1 

ISTOP=1    =>  STOP  FPUV  UPON  RETURN  TO  MAIN  PROG 

NSR        =>  NUMBER  TIME  STEPS  IN  REVOLUTION 

NR  *>  RADIUS  NUMBER 

XFP.YBAR.ZBAR.RFP.TZFP 

=>  X.Y.Z. RADIAL. THETA  FIELD  PT  COORD 
THETA(1-5)=>  ANGLES  5  TIME  STEPS  BEFORE  ZERO 

BLAOE  ANGLE 
THETA(6  -  NSR+5)  =>  ANGLES  IN  ONE  BLADE  REV. 
THETA(NSR+6  -  NSR+20)  =>  ANGLES  15  TIME  STEPS 

AFTER  ZERO  BLAO  ANGLE. 
VX(K)       =>  AXIAL  VELOCIT?  AT  ANGLE  THETA(K) 
UA(KK.M)   =>  SMOOTH  AXIAL  VELOCITY  AT  RADIUS  M. 

ANGLE  KK 


***  COMMON  TO  CONTAIN  GEOMETRIC  PARAMETERS  , 

COMMON  /  GEOCOM  /  MM, MN, NN, NW, NTORL .X( 12 . 10) . Y( 12 . 10) , Z( 12 . 10) . 
&    XTIP( 10,  11).YTIP( 10. 11) ,ZTIP(  10.  1 1 ).XH( 101 ,7),YH( 101,7) , 
&    ZH( 101.7) ,XW( 21 , 10).YW(21 . 10) ,ZW(21 . 10) . RH 

«**  COMMON  TO  CONTAIN  SINGULARITY  STRENGTHS 

COMMON  /  SINCOM  /  GT( 10. 10. 60) . G8( 10. 9 . 60) . GTW( 10) . SB ( 10. 9 ) . GTV 
&   .GAM(9,60) ,GMEAN(9),DKFW(20) .HBWX . HBWY . HBWZ 

•*-  COMMON  TO  CONTAIN  FIELD  POINT  DATA.  N8LADE .  IDENT.  NSR.  AND  VSR 

COMMON  /  FPCOM  /  N8LADE . NSR . XFP , RFP . TZFP , YBAR . ZBAR . 
&    DBLAOE.DTFP. I DENT ( 18) ,TFP(60) .VSR. I  PLOT. 01 AM 

»•*  COMMON  TO  CONTAIN  VELOCITY  DATA 

COMMON  /  VELCOM  /  VIX(60) , VIT( 60) , VIR( 60) .U( 60, 3) . NSAMP 

***  COMMON  TO  CONTAIN  HARMONICS  COEFFICIENTS 


UFP07460 
UFP07470 
UFP07480 
UFP07490 
UFP075O0 
UFP07510 
UFP07520 
UFP07530 
UFP07540 
UFP07550 
UFP07560 
UFPO7570 
UFP07580 
UFP07590 
UFP07600 
UFP07610 
UFP07620 
UFP07630 
UFP07640 
UFP07650 
UFP07660 
UFP07670 
UFP07680 
UFP07690 
UFP0770O 
UFP077 10 
UFP07720 
UFP07730 
UFP07740 
UFP07750 
UFP07760 
UFP07770 
UFP07780 
UFP07790 
UFP07800 
UFP07810 
UFP07820 
UFP07830 
UFP07840 
UFP07350 
UFP07860 
UFP07870 
UFP07880 
UFP07890 
UFP07900 
UFP07910 
UFP07920 
UFP07930 
UFP07940 
UFP07950 
UFP07960 
UFP07970 
UFP07980 
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c 

c  •  *  * 

c 

102 

C 

c  *** 

c 
c 

c  »  *  • 


c 

c  «*- 
c  *«* 

c 

200 

c 

c  **- 

c 

300 
C 

c  *** 

c 


c 

c  * ** 

c 


100 

c 

c  »»* 

c 


COMMON  /  HARCOM  /  A( 15 , 3) ,B( 15 . 3 ) , AMP( 15 , 3 ) , PH( 15 , 3 ) ,NH 

DIMENSION  THETA(aO) . VX( 80) , UA ( 60, 9 ) , RWK(9) 

N=1 

NR=1 

CALCULATE  VELOCITIES  AT  FIVE  BLAOE  GRID  RADII 

DO  300  M  =  3.9.2 

NR=NR+1 

CALCULATE  VELOCITIES  AT  NSR  ANGLES  FOR  EACH  RADIUS 

DO  200  KK=1 .NSR 

FIELD  POINT  CENTERED  ON  BLADE  GRID 

CALL  FPGRID(X. Y.Z.N, M.KK, 12) 

IF  (M.E0.3)  X1=XFP 

KK15=KK+15 

THETA(KK15)=TZFP 

RWK(NR)=RFP 

DTFP=360. /FLOAT (NSR) 

CALCULATE  AXIAL  VELOCITY  FOR  FIELD  POINT  AT  RADIUS  RWK(NR)  AND 
ANGLE  THETA(KK15).   STORE,  THIS  VELOCITY  IN  VX . 

CALL  AMPCAL(VX( 16) . KK ) 

NSR15=NSR+15 

STORE  SMOOTHED  VELOCITIES  IN  NR-TH  COLUMN  OF  UA 

CALL  SMOOTH(NSR. NR.VX.THETA.UA) 

CALCULATE  NSR  VELOCITIES  AT  PROPELLER  HUB 

NRHUB=1 
XFP=1  .5*X1 
RWK(NRHUB)=RH- 

STEP  THROUGH  NSR  ANGLES  AT  HUB 

DO  10O  KK=1 .NSR 

KK15=KK+15 

THETA(KK15)=360. «FL0AT(KK-1 )/FLOAT(NSR) 

TZFP=THETA(KK15) 

CALL  AMPCAL(VX( 16) , KK ) 

NSR15=NSR+15 

NOW  STORE  SMOOTHED  HUB  VELOCITIES  IN  FIRST  COLUMN  OF  UA 


UFP07990 
UFP080CO 
UFP08010 
UFP08020 
UFP08030 
UFP08040 
UFPCS050 
UFP08060 
UFP08070 
UFP08080 
UFP08090 
UFP08  100 
UFP08 1  10 
UFP08  120 
UFP08130 
UFP08 140 
UFP08  150 
UFP08 160 
UFP08170 
UFP08180 
UFP08  190 
UFP08200 
UFPO8210 
UFP08220 
UFP08230 
UFP08240 
UFP08250 
UFP08260 
UFP08270 
UFP08280 
UFP08290 
UFP08300 
UFP08310 
UFP08320 
UFP08330 
UFP08340 
UFP08350 
UFP08360 
UFP08370 
UFP08380 
UFP08390 
UFP08400 
UFP08410 
UFP08420 
UFP08430 
UFP08440 
UFP08450 
UFP08460 
UFP08470 
UFP08480 
UFP08490 
UFP0850O 
UFP08510 
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c 

c  •»« 

c  *** 

c 

c  &&& 

c 
c 
c 

c  *** 


c 

c  *«■ 

c 


400 

c 

c  **• 

c 

500 

c 

c  * *« 

c 

901 
902 
101 


C 
c  *» 

c 


CALL  SM00TH(NSR . NRHU8 , VX . THETA , UA ) 

NOW  CALCULATE  VELOCITIES  AT  A  NEXT  RADII  OUTSIDE  THE  PROPELLER 
DISK.   USE  XFP  OF  PROPELLER  GRID  CASE  OF  N= 1 .  M  =  9. 

DIAGNOSTICS  &&&S&&&&&&&&&&&&S 

GOTO  101 

CALCULATE  NSR  VELOCITIES  AT  1.1.1.4.  AND  1.7  PROP  RADII 

DO  500  KKK=1 .7.3 
RFP=1 . 0+FLOAT(KKK)/lO. 
NR=NR+1 
RWK(NR)=RFP 

CALCULATE  NSR  VELOCITIES  AT  THESE  RADII 

DO  400  KK=1 .NSR 

KK15=KK+15 

THETA(KK15)=36O.*FL0AT(KK-1 )/FLOAT(NSR) 

TZFP=THETA(KK15) 

CALL  AMPCAL(VX( 16) ,  KK  ) 

NOW  STORE  SMOOTHED  VELOCITIES  IN  NR-TH  COLUMN  OF  UA 

CALL  SMOOTH(NSR. NR.VX.THETA.UA) 

WRITE(OUTPUT  IN  FORMAT  FOR  UNSTEADY  PIEWAKE 

F0RMAT(8F10.5) 

F0RMAT( 1615) 

WRITE(8.902)  NR 

WRITE(8.901)  (RWK(M) .M=1 . NR ) 

WRITE (8. 901 )  ( (UA(KK.M).KK=1 .NSR) . M=1 . NR ) 

*  STOP  FPUV  UPON  RETURN  TO  MAIN  PROGRAM 

IST0P=1 
RETURN       / 
END 


UFP08520 
UFP08530 
UFP08540 
UFP08550 
UFP08560 
UFP08570 
UFP08580 
UFP08590 
UFP086OO 
UFP08610 
UFP08620 
UFP08630 
UFP08640 
UFP08650 
UFP08660 
UFP08680 
UFP08690 
UFP0870O 
UFP087  10 
UFP08720 
UFP08730 
UFP08740 
UFP08750 
UFP08760 
UFP08770 
UFP08780 
UFP08790 
UFP08800 
UFP08810 
UFP08S20 
UFP08830 
UFP08840 
UFP08350 
UFP08860 
UFP08870 
UFP08880 
UFP0889O 
UFP08900 
UFP08910 
UFP08920 
UFP08930 
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SUBROUTINE  AMPCAL( VX , KK ) 


«******#***J 


*********** 


I  *****  * 


SUBROUTINE  AMPCAL  WILL  TAKE  FIELD  POINT  GEOMETRY 
CALCULATED  IN  PIEDAT  ANO  CALCULATE  THE  INDUCED 
TIME-AVERAGED  AXIAL  VELOCITY.   IT  WILL  STORE  THAT 
VELOCITY  IN  VX(KK) . 

NSR        =>  NUMBER  OF  ANGLES  PER  REVOLUTION 
K  =>  ANGLE  INCREMENTS  +  1 

VX(KK)     »>  WHERE  IME-AVERAGED  AXIAL  VEL  STORED 
AMP(1,1)   =>  ZEROETH-HARMONIC  AXIOAL  VELOCITY 
COMPUTED  BY  FOUCAL 


C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 

c 

C  «*«  COMMON  TO  CONTAIN  GEOMETRIC  PARAMETERS 
C 

COMMON  /  GEOCOM  /  MM . MN , NN , NW . NTORL . X( 12 .  10) . Y (  1 2 .  10 ) . Z (  1  2 .  10) . 
&    XTIP(  10.  1 1 ) ,YTIP( 10.  1  1  ).ZTIP(  10, 11 ) .XH( 101 .7) . YH( 101 ,7) , 
&    ZH( 101 ,7).XW(21 . 10). YW(21. 10) ,ZW(21 . 10) ,RH 
C 

C  •**  COMMON  TO  CONTAIN  SINGULARITY  STRENGTHS 
C 

COMMON  /  SINCOM  /  GT ( 10 . 10 . 60) . GB ( 10 . 9 . 60 ) . GTW( 10) , SB ( 10 . 9 ) . GTV 
&   , GAM (9, 60) ,GMEAN(9) . DKFW(20) . HBWX . HBWY . HBWZ 
C 

C  »«»  COMMON  TO  CONTAIN  FIELD  POINT  DATA.  NBLADE.  IDENT.  NSR.  AND  VSR 
C 

COMMON  /  FPCOM  /  NBLADE . NSR . XFP . RFP . TZFP . YBAR . ZBAR . 
&    DBLADE,DTFP,IDENT(  18) ,TFP(60) .VSR. I  PLOT. D I AM 
C 

C  ***  COMMON  TO  CONTAIN  VELOCITY  DATA 
C 

COMMON  /  VELCOM  /  VIX( 60) . VIT( 60) . VIR( 60) . U( 60. 3 ) . NSAMP 
C 

C  »•*  COMMON  TO  CONTAIN  HARMONICS  COEFFICIENTS 
C 

COMMON  /  HARCOM  /  A ( 1 5 . 3 ) . B ( 1 5 , 3 ) . AMP ( 1 5 . 3 ) . PH( 1 5 , 3 ) . NH 
C 

DIMENSION  VX(  1  ) 
C 

C  »•*  RUN  THROUGH  AND  FPSTEP  FOR  EACH  TIME  STEP 
C 

DO  100  K=1 ,NSR 
C 

C  *•»  FPSTEP  CALCULATES  VELOCITY  INDUCED  AT  A  FIELD  POINT  BY  ONE 
C  »**  BLADE  AT  TIME  STEP  K 
C 

100   CALL  FPSTEP(K) 
C 


UFP08940 
UFP08950 
UFP08960 
UFP08970 
UFP08980 
UFP08990 
UFP09OO0 
UFP09010 
UFP09020 
UFP09030 
UFP09040 
UFP090S0 
UFP09060 
UFP09070 
UFP09080 
UFP09090 
UFP0910O 
UFP091 10 
UFP09120 
UFP09  130 
UFP09140 
UFP09150 
UFP09  160 
UFP09170 
UFP09  180 
UFP09  190 
UFP0920O 
UFP092  10 
UFP09220 
UFP09230 
UFP09240 
UFP09250 
UFP09260 
UFP09270 
UFP09280 
DFP09290 
UFP09300 
UFP09310 
UFP09320 
UFP09330 
UFP09340 
UFP09350 
UFP09360 
UFP09370 
UFP09380 
UFP09390 
UFP0940O 
UFP09410 
UFP09420 
UFP09430 
UFP09440 
UFP09450 
UFP09460 
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C  **»  SUMVEL  SUMS  FIELD  POINT  VELOCITIES  INDUCED  BY  EACH  BLADE  UFP09470 

C  UFP09480 

CALL  SUMVEL  UFP09490 

C  UFP09500 

C  -**  FOUCAL  GENERATES  HARMONIC  COEFFICIENTS  FOR  INDUCED  FIELD  POINT     UFP09510 

C  •*•  VELOCITY  AS  PROPELLER  ROTATES.  UFP09520 

C  UFP09530 

CALL  FOUCAL  UFP09540 

C  UFP09550 

C  •••  STORE  TIME-AVERAGED  ( ZERO -HARMONI C )  VELOCITY  IN  VX  UFP09560 

C  UFP09570 

VX(KK)=AMP( 1 , 1 )  UFP09580 

RETURN  UFP09590 

END  '                       UFP09600 


/ 
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SUB ROUT  I NE  SMOOTH ( NSR .  NR , VX . THETA , UA ) 


SUBROUTINE  SMOOTH  REAOS  INDUCED  VELOCITIES  AT  ONE 
RADIUS  FROM  PIEDAT.   IT  CALCULATES  BY 
INTERPOLATION  A  SERIES  OF  NSR  EVENLY  SPACED 
AXIAL  INOUCED  VELOCITIES  SUITABLE  FOR  INPUT  INTO 
UNSTEADY  PIEWAKE. 


NSR         =>  NUMBER  OF  ANGLES  PER 
NR  =>  RADIUS  NUMBER  1=HUB. 

VX(1  -  5)  =>  RAW  VELOCITIES  5  STEP 

BLADE  ANGLE. 
VX(6  -  NSR+5)  =>  RAW  .VELOCITIES  IN 
VX(NSR+6  -  NSR+20)  =>  RAW  VELOCITI 

TIME  STEPS  AFTER  ONE 
THETA(K)  =>  ANGLE  OF  VELOCITY  VX( 
UA(K.M)  =>  SMOOTH  VEL  AT  RADIUS 
TH(K)  =>  EVENLY  SPACED  ANGLES 
NIN  =>  NUMBER  OF  ANGLES  FED 
KPOS  =>  K  OF  FIRST  POS.  ANGLE 
KSTART.KSTOP  *>  FIRST  AND  LAST  IND 

AND  VX(K)  SENT  TO  UGL 
AVX        =>  SPLINE  CURVE  COEFFICI 


REVOLUTION 

ETC. 

S  BEFORE  ZERO 

ONE  REVOLUTION 
ES  15  TIME 
REVOLUTION 
K) 

M.  ANGLE  K 
OF  UA(K.M) 
TO  UGLYDK 

IN  THETA(K) 
EX  OF  THETA(K) 
YDK 
ENTS  OF  VX 


C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
•c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

DIMENSION  VX(80) ,UA(60,9)  .THETA (80) .TH(60) , AVX (250) 
C 

C  •«*  COPY  LAST  15  STEPS  IN  REV  INTO  FIRST  15  ENTRIES  IN  VX  AND  THETA 
C 

DO  100  K=1 , 15 

KNSR=K+NSR 

THETA (K)=THETA(KNSR) -360. 
100   VX(K)=VX(KNSR) 
C 

C  •»*  COPY  FIRST  5  TIME  STEPS  IN  REV  INTO  5  ENTRIES  AFTER  REVOLUTION 
C 

DO  200  K=16,20 

KNSR=K+NSR 

THETA(KNSR)=THETA(K)  +  360. 

VX(KNSR)=VX(K) 

SET  UP  EVENLY  SPACED  ANGLES 

DO  300  K=1 .NSR 

TH(K)=360. -  FLOAT (K- 1 ) /FLOAT (NSR) 

FORMAT(8F10.5) 

FORMAT( 1615) 

FIND  FIRST  POSITIVE  ANGLE.  SAVE  INDEX 


200 
C 
C  •»* 

c 

300 
901 
902 

C 

c  »»* 

c 


UFP096  10 
UFP09620 
UFP09630 
UFP09640 
UFP09650 
UFP09660 
UFP09670 
UFP09630 
UFP09690 
UFP09700 
UFP097  10 
UFP09720 
UFP09730 
UFP09740 
UFP09750 
UFP09760 
UFP09770 
UFP09780 
UFP09790 
UFP0980O 
UFP09810 
UFP09820 
UFP09830 
UFP09840 
UFP09850 
UFP09860 
UFP09870 
UFP09880 
UFP09890 
UFP09900 
UFP09910 
UFP09920 
UFP09930 
UFP09940 
UFP09950 
UFP09960 
UFP09970 
UFP09980 
UFP09990 
UFP100O0 
UFP10010 
UFP10020 
UFP10030 
UFP10040 
UFP10050 
UFP10060 
UFP10070 
UFP10080 
UFP10090 
UFP10100 
UFP101 10 
UFP10120 
UFP10130 
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DO  40O  K=1 .NSR 

KPOS=K 

IF  (THETA(K) .GT.O. )  GOTO  500 
400   CONTINUE 
500   KSTART=KPOS- 1 

NSR2=NSR/2 

KST0P=KP0S+NSR2+ 1 

NIN=KST0P-KSTART+1 

:  «*•  FINE  SPLINE  COEFFICIENTS  AND  SMOOTH  VEL  FOR  FIRST  180  DEGREES 

CALL  UGLYDK(NIN, 1 . 1 . THETA (KSTART ) , VX ( K START ) . 0 . ,0. , AVX) 
CALL  EVALDK(NIN.NSR2, THETA (KSTART) . TH( 1 ) , UA ( 1 , NR ) . AVX ) 


C 

c  *** 

C 


FIND  SPLINE  COEFFICIENTS  AND  SMOOTH  VEL  FOR  LAST  180  DEGREES 

KSTART=KST0P-2 

KSTOP=KSTART+NIN- 1 

NSR2P=NSR2+1 

CALL  UGLYDK(NIN. 1 . 1 . THETA ( KSTART ) ,VX( KSTART) ,0. .0. .AVX) 

CALL  EVALDK(NIN.NSR2, THETA ( KSTART ).TH(NSR2P) . UA (NSR2P . NR ) ,AVX) 

RETURN 

END 


UFP10140 
UFP10150 
UFP10160 
UFP10170 
UFP10180 
UFP10190 
UFP10200 
UFP102 10 
UFP10220 
UFP10230 
UFP10240 
UFP10250 
UFP10260 
UFP10270 
UFP10280 
UFP10290 
UFP1030O 
UFP10310 
UFP10320 
UFP10330 
UFP 10340 
UFP10350 
UFP 10360 


/, 
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**  APPENDIX  XI   ** 


PIEVAXE  Program  Listing 


j? 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


THIS  PROGRAM  COMPUTES  EFFECTIVE  WAKES  BY  DNV  METHOD 


THE  ORIGINAL  PIEWAKE  WAS  WRITTEN  BY  R.  VAN  HOUTEN 
IN  ABOUT  1981.   IT  IS  BASED  ON  THE  DET  NORSKE 
VERITAS  METHOD  OF  COMPUTING  EFFECTIVE  WAKE.   THIS 
PROGRAM  HAS  BEEN  MODIFIEC  BY  R.  JAMISON  IN  APRIL. 
1982.  TO  ACCOMODATE  UNSTEADY  PROPELLER  -  I NDUCED 
AXIAL  VELOCITIES  FROM  UFPV.   BASICALLY.  BOTH 
PROGRAMS  DIVIDE  THE  WAKE  INTO  A  NUMBER  OF  PIE- 
SHAPED  SEGMENTS.   EACH  SEGMENT  IS  ASSUMED  TO 
CONTRACT  IN  AN  AXISYMMETRIC  MANNER  AS  GIVEN  BY 
T.  HUANG. 


NP 

NX 

RX(M) 

VO(K.M) 

VT(K.M) 

VR(K.M) 

VTE. VRE 

UASEG(M) 

UX(M) 


AUX. AUASEG=> 


=  > 

=  > 
=  > 
=  > 
=  > 
=  > 
=  > 


NUMBER  OF  PIE-SHAPED  SEGMENTS 

NUMBER  OF  NOMINAL  WAKE  RADII 

NOMINAL  WAKE  RADIUS  M  (UP  TO  11) 

NOMINAL  AXIAL  VELOCITY.  SEG  K.  RAD  M 

NOM  TANG  VEL,  SEG  K.  RAD  M 

NOM  RAD  VEL.  SEG  K.  RAD  M 

VT  &  VR  INTERPOLATED  AT  EFF  WAKE  RAD 

UA  IN  A  PIE-SEGMENT.   ALSO  USED  TO 

STORE  A  PIE  -SEGMENT  OF  VT 

VO  IN  A  PIE-SEGMENT.   ALSO  USED  TO 

STORE  A  PIE-SEGMENT  OF  VR 

SPLINE  COEFF.  FOR  VR . VT .  RESP. 


RE(M) 

THETA(K) 

NA 

RA(M) 

UA(K.M) 

VE(K.M) 

VL 

VLAVN 

VLAVE 

BLOCK 


=  > 
=  > 
=  > 
=  > 
=  > 
=  > 

=  > 
=  > 


RADII 
ANGLE 


OF  EFFECTIVE 
IN  DEGREES  OF 


NUMBER  OF  PROP-INO, 
PROP-IND  AXIAL  VEL 
PROP  IND  AXIAL  VEL, 
EFFECTIVE  WAKE  VEL, 
AVERAGE  VELOCITY 
VOLUMETRIC  AVERAGE 
VOLUMETRIC  AVERAGE 
EFFECTIVE  BLOCKAGE 


WAKE 

VELOCITIES 

PIE 

SEG  K 

AXIAL 

VEL  RADI 

RADIUS 

M 

SEG 

K 

RAD  M 

SEG 

K 

RAO  M 

NOMINAL  VELOCITY 
EFFECTIVE  VELOCIT 


DIMENSION  VO(60. 1  1  ) . RX( 1 1 ) . UX ( 60 ) . RA (  1 1 ) , UA ( 60 ,  1  1  )  . UASEG( 80) 
DIMENSION  RE(  1 1) .  VR( 60. 1 1 ) , VT( 60, 1  1  )  , AUASEG( 44 ) , AUX (44) ,VTE(  1 1 ) 
DIMENSION  THETA(60) . VE(60. 11 ) ,IDENT( 18) .UERX(60) ,VRE( 1 1 ) 


1  10 

103 

1  13 

1  14 

115 

111 

210 

20O 

c 

c 

*  *  * 

F5.3) 
''  .F5.3) 


COMMON  RX.UX.RA.UASEG.UERX.RE 

F0RMAT(3F 10.5) 

FORMAT(/'  AVERAGE  VELOCITY  * 

FORMAT(/'  VOLUMETRIC  AVERAGE 

FORMAT(/'  VOLUMETRIC  AVERAGE 

FORMAT(/'  EFFECTIVE  BLOCKAGE 

FORMAT( 18A4) 

FORMAT* 1615) 

FORMAT(8F10.5) 


REAO  IN  NOMINAL  WAKE  OATA.   SAVE  RADIAL  AND  TANGENTIAL 


.F5.3) 

NOMINAL  VELOCITY 
EFFECTIVE  VELOCITY 
«'.F7.3.'  PERCENT') 


PIE00010 
PIE0O020 
PIE0O030 
PIE00040 
PIE0O050 
PIEOO060 
PIEOOO70 
PIE00O80 
PIEOO090 
PIE0O10O 
PIE0O1 10 
PIE00120 
PIE00130 
PIE00140 
PIE00150 
PIE0O160 
PIE0O170 
PIE0O180 
PIE00190 
PIE002O0 
PIE0O210 
PIE00220 
PIE00230 
PIE0O240 
PIE0O250 
PIE00260 
PIE0O270 
PIE00280 
PIE0O290 
PIE00300 
PIE003  10 
PIE0O320 
PIE0O330 
PIEOO340 
PIE00350 
PIE0O360 
PIECO370 
PIE00380 
PIE00390 
PIE0O40O 
PIECO410 
PIE00420 
PIE0O430 
PIE0O440 
PIE0O4  50 
PIE0O460 
PIE0O4  70 
PIEOO480 
PIE0O4  90 
PIE00500 
PIE0O510 
PIE00520 
PIEOOS30 
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VELOCITIES.   THEY  ARE  NOT  NEEDED  FOR  EFFECTIVE  WAKE 
CALCULATIONS  BUT  WILL  BE  ECHOED  FOR  WKPROC. 

READOO.111)  (  IDENT(  I  )  .  1  =  1  .  18) 
READ( 10.210)  NX 
READ( 10.210)  NP 
READ(  10. 200) (RX(M) ,M= 1 .NX) 
DO  20  M=1 .NX 

REAO(  10. 200) (  (THETA(K) . VO ( K  ,  M )  . VT ( K . M ) , VR ( K . M ) ) , K= 1 . NP ) 
20  CONTINUE  •  . 

"  READ  IN  PROPELLER-INDUCED  VELOCITIES  FROM  EITHER  FPUV 
OR  A  SEPARATE  DATA  FILE 


900   F0RMAT(/'  DO  YOU  WISH  TO  USE  STEADY  OR  UNSTEADY  VELOCITY  PROFILE? 
&   /'  ENTER  1  FOR  STEADY,  2  FOR  UNSTEADY') 
8   WRITE(6.90O) 

READ(5.»)  ISTED 

IF  (ISTED. GT. 2. OR. ISTED. LT. 1 )  GOTO  8 
READ(8.210)  NA 
READ(8.20O)  ( RA ( I ) . I = 1 , NA ) 
RH  =  RA(  1  ) 
RPR0P=1  . 

IF(ISTED.EQ. 1 )  GOTO  22 

READ(8.20O)  ( ( UA ( KK . M ) , KK= 1 . NP ) . M= 1 . NA ) 
GOTO  92 
22   REA0(8,2OO)  ( UA( 1 . M) , M= 1 . NA ) 

DO  21  KK=2.NP 
*    DO  21  M=1 ,NA 
21   UA(KK,M)=UA( 1 ,M) 
C 

C  •**  EXTRAPOLATE  NOMINAL  WAKE  DATA  TO'  HUB 
C 

92   WRITE(6.10O) 

100  FORMAT(/'  DO  YOU  WISH  TO  EXTRAPOLATE  VELOCITY  DATA  TO  THE  HUB?'/ 
1   '  ENTER  2  FOR  NO.  1  FOR  LINEAR,  O  FOR  CONSTANT') 
READ(5. * )IXTRAP 
NX1=NX 

IF(IXTRAP.GT. 1 )G0  TO  28 
C 

C  ***  LINEAR  EXTRAPOLATION  OF  NOMINAL  WAKE  DATA  TO  HUB. 
C 

NX=NX+1 
DO  25  K=1 .NP 
DO  23  M=1 ,NX1 
MR=NX-M 
MRP1=MR+1 

VT(K.MRP1 )=VT(K.MR) 

VR(K.MRP1 )*VR(K,MR) 

23  V0(K.MRP1 )=VO(K,MR) 

25  VO(K. 1 )=V0(K.2)+IXTRAP*(V0(K.2)-V0(K,3))/(RX( 1 ) -RX( 2 ) ) * ( RH-RX ( 1 ) ) 
DO  26  M=1 ,NX1 
MR=NX-M 


PIE00S40 
PIE0O55O 
PIE00560 
PIE00570 
PIE0O580 
PIE00590 
PIE00600 
PIE006  10 
PIE00620 
PIE0O630 
PIE00640 
PIE00650 
PIE0O660 
PIE00670 
'PIE006S0 
PIE0O690 
PIEO070O 
PIECO710 
PIEOO720 
PIE00730 
PIE00740 
PIE00750 
PIE00760 
PIE00770 
PIE0O780 
PIE0O790 
PIE0080O 
PIE0O810 
PIE0O820 
PIE0O830 
PIE00840 
PIE00850 
PIE0O860 
PIE0O870 
PIE0O880 
PIE00890 
PIE009CO 
PIE00910 
PIE00920 
PIE0O930 
PIE00940 
PIE00950 
PIE00960 
PIECO970 
P.IE0O980 
PIE00990 
PIE010OO 
PIE01010 
PIE01020 
PIE01030 
PIE01040 
PIE01O5O 
PIE01060 
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MRP1=MR+1  PIE01070 

26  RX(MRP1 )=RX(MR)  PIE01Q80 

RX(1)=RH  PIE01090 

28  VL=0.  PIE01 100 

VLAVN=0.  PIE01110 

VLAVE=0.  PIE01120 

C  PIE01 130 

C  •**  CYCLE  THROUGH  EACH  PIE  SEGMENT  PIE01140 

C  PIE01150 

DO  50  K=1 .NP  PIE01 160 

DO  30  M=1 .NX  PIE01 170 

UASEG(M)=UA(K.M)  PIE01180 

30  UX(M)=VO(K.M)  PIE01190 

C  PIE01200 

C  -««  WKMOD  PERFORMS  THE  HUANG  MODIFICATION  ON  ONE  PIE-SEGMENT  PIE01210 

C  PIE01220 

CALL  WKMOD (NX. NX. NA, 1 .RPROP, VOL. VOLAVN, VOLAVE)  PIE01230 

C  PIE01240 

C  *«*  INTERPOLATE  VT  AND  VR  AT  EFFECTIVE  WAKE  RADII  PIE01250 

C  PIE01260 

DO  60  M= 1 .NX  PIE01270 

UASEG(M)=VT(K.M)  PIE01280 

60   UX(M)=VR(K,M)  PIE01290 

CALL  UGLYDK(NX.  1 .  1  .RX.UASEG.O.  .0.  . AUASEG)  PIE01300 

CALL  UGLYDK(NX. 1 . 1.RX.UX.0. .0. .AUX)  PIE01310 

CALL  EVALDK(NX. NX. RX. RE. VTE. AUASEG)  PIE01320 

CALL  EVALDK(NX,NX.RX,RE. VRE, AUX)  PIE01330 

DO  70  M=1,NX  PIE01340 

VT(K.M)=VTE(M)  PIE01350 

70   VR(K.M)=VRE(M)  PIE01360 

C  PIE01370 

C  •*-  ACCUMULATE  VOLUMETRIC  AVERAGE  DATA  PIE01380 

C  PIE01390 

VL=VL+VOL                                -  PIE01400 

VLAVN=VLAVN+VOLAVN  PIE014  10 

VLAVE=VLAVE+VOLAVE  PIE01420 

DO  40  M=1 .NX  PI EO 1430 

40  VE(K.M)=UERX(M)         '  PIE01440 

50  CONTINUE  PIE01450 

C  PIE01460 

C  ••-  COMPUTE  VOLUMETRIC  AVERAGES  PIE01470 

C                                                      .  PIE01480 

VL=VL/NP  PIE01490 

VLAVN=VLAVN/NP  PIE015OO 

VLAVE=VLAVE/NP  PIE01510 

SLOCK =( 1 . -VLAVE/VLAVN)* 100.  PIE01520 

C  PIE01530 

C  ...  WRITE  EFFECTIVE  WAKE  DATA  IN  FILE  11  TO  BE  PROCESSED  BY  WKPROC     PIE01540 

C  PIE01550 

WRITE( 11 .  111)  (IDENT(I).I-I.  18)  PIE01560 

WRITEC 1 1 .210)  NX  PIE01570 

WRITE( 1 1 .210)  (NP.M=1.NX)  PIE01580 

WRITE( 1 1 ,200)(RE(M) .M=1 .NX)  PIE01590 
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DO  120  M=1.NX  PIE016OO 

WRITE(  1  1  .200)( (THETA(K)  ,VE(K,M).VT(K.M).VR(K,M)),K=1,NP)  PI EO 16 10 

120  CONTINUE  PIE01620 

C  P1E01630 

C  •»•  WRITE  AVERAGES  TO  TERMINAL                                           PIE01640 

C  PIE01650 

WRITE(6. 103)VL  PIE01660 

WRITE(6, 1  13)VLAVN  PIE01670 

WRITE(6.  1  14)VLAVE  PIE01680 

WRITE(6.  1  15)BL0CK  PIE01690 

STOP  PIE0170O 

END  PIE01710 
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SUBROUTINE  WKMOO ( N , NX . NA , I OW , RPROP , VOL . VOLAVN , VOL AVE ) 


C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


c 

c  *«* 

c 


C   10 
10 

c  , 

Q     »  *  * 

c 


SUBROUTINE  WKMOO  CALCULATES  THE  THOMAS  HUANG 
CONTRACTION  OF  THE  NOMINAL  AND  INDUCED  WAKES 
PASSED  TO  IT  THROUGH  COMMON.   IT  ALSO  COMPUTES 
VOLUMETRIC  AVERAGES. 


NX 

=  > 

RH 

=  > 

RT 

=  > 

UX.UXR 

=  > 

UE .UEXR 

=  > 

UP 

=  > 

UA.UAR 

=  > 

IOW  =  1 

=  > 

ITER 

=  > 

ITERAT 

=  > 

B.C.D.F 

=  > 

AUA. AUX. 

AUE  = 

RPROP 

=  > 

RE(M) 

=  > 

VOL 

=  > 

VLAVN 

=  > 

VLAVE 

=  > 

NUMBER  OF  NOM 
RADIUS  OF  HUB 
RADIUS  OF  OUT 
NOMINAL  WAKE 
EFFECTIVE  WAK 
APPARENT  WAKE 
INDUCED  VELOC 
NO  TUNNEL  COR 
OPEN  WATER  IT 
HUANG  CORRECT 
HUANG  CORRECT 
>  SPLINE  COEFF 
AND  EFFECTIVE 
RADIUS  OF  PRO 
SPECIFIED  EFF 
AVERAGE  VELOC 
VOLUMETRIC  AV 
VOLUMETRIC  AV 


INAL  RADII 

ERMOST  NOMINAL  VELOCITY 

VELOCITIES 

E  VELOCITIES 

=  EFF.  WAKE  +  INO.  WAK 
ITIES 
RECTIONS 
ERATION  INDEX 
ION  FACTOR  INDEX 
ION  FACTORS.  SEE  PAPER 
.  FOR  INDUCED.  NOMINAL. 

VELOCITY  FIELDS.  RESP 
PELLER 

ECTIVE  WAKE  RADII 
ITY 

ERAGE  NOMINAL  VELOCITY 
ERAGE  EFFECTIVE  VELOCITY 


DIMENSION  RX( 1 1).R(  1  1  )  ,RA( 1 1 ) .UX(60) . UXR( 60 ) . UA( 60 ) .UP (60) 
DIMENSION  RE( 1 1 ).UAR(60) ,RP( 11 ) 
DIMENSION  AUX(240) ,AUA(240) 
DIMENSION  AUE(240) ,UERX(60) 
COMMON  RX.UX.RA,UA„UERX.RE 

*-*  USE  R(I)  FOR  NOMINAL  RADII 


N  =  NX 
N1=N-1 
RH=RX( 1 ) 
RT=RX(NX) 

INTERPOLATE  TO  FIND  NOMINAL  VELOCITIES  AT  NOMINAL  RADII 

CALL  UGLYDK(NX, 1 , 1 .RX.UX.O. .0. .AUX) 

DO  10  1  =  1  ,N 

R(I)=(RT-RH)/(N-1 )«(I-1 )+RH 

R(I)=RX(I) 

CALL  EVALDK(NX.N.RX.R.UXR.AUX) 

INTERPOLATE  TO  FIND  INDUCED  VELOCITIES  AT  NOMINAL  RADII 

CALL  UGLYOK(NA. 1 . 1 . RA . UA . 0 . . 0 . .AUA) 


UE(60) 


PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 
PI 


E01720 
E01730 
E01740 
E01750 
E01760 
E01770 
E01780 
E01790 
E0180O 
E01810 
E01820 
E01830 
E01840 
E01850 
E01860 
E01870 
E01880 
E01890 
E01900 
E01910 
E01920 
E01930 
E01940 
E01950 
E01960 
E01970 
E01980 
E01990 
E02OO0 
E02010 
E02020 
E02030 
E02040 
E02050 
E02060 
E02070 
E02080 
E02090 
E02  1CO 
E021 10 
E02  120 
E02130 
E02140 
E02150 
E02160 
E02  170 
E02 180 
E02  190 
E02200 
E02210 
E02220 
E02230 
E02240 


-  94  - 


CALL  EVALDK(NA,N.RA ,R,UAR, AUA )  PIE02250 

901   F0RMAT(8F 10. 5)  PIE02260 

C  PIE02270 

C  •«*  OUTSIDE  PROPELLER  DISK.  NOMINAL  WAKE  =  EFFECTIVE  WAKE  PIE02280 

C  PIE02290 

IUPDAT=0  PIE0230O 

ITER=1  PIE02310 

ITERAT=1  PIE02320 

UE(N)=UXR(N)  PIE02330 

20  UP(N)=UE(N)+UAR(N)  PIE02340 

C                                             n  PIE02350 

C  •*•  ESTIMATE  EFF.  WAKE  MARCHING  FROM  OUTSIDE  TO  HUB  PIE02360 

C  PIE02370 

DO  30  1=1, N1  PIE02380 

K=N-I  PIE02390 

UE(K)=S0RT( (UE(K+1  )  +  (UAR(K+1 )+UAR(K) )/2.  ) • • 2  +  UXR ( K ) « » 2 -UXR ( K+ 1 )»»2PIE02400 

1 )-(UAR(K+ 1 )+UAR(K) )/2.  PIE02410 

UP(K)=UE(K)+UAR(K)  PIE02420 

30  CONTINUE  PIE02430 

C  PIE02440 

C  *"*  COMPUTE  HUANG'S  CORRECTION  FACTORS  PIE02450 

C  PIE02460 

RP(1)=RH  PIE02470 

VOL=0.  PIE02430 

DO  40  K=1 ,N1  PIE02490 

B=2«UP(K+1 )+UP(K)  PIE02500 

C=-RP(K)«(UP(K+1 )-UP(K))  PIE02510 

F  =  (R(K+1 )**2-R(K)**2)«(2«UXR(K+1 )+UXR(K) ) -R ( K  )  « ( UXR ( K+ 1 )-UXR(K) )*(PIE02  520 

1R(K+1 )-R(K))  PIE02530 

V0L=V0L+F  PIE02540 

0=-RP(K)**2*(UP(K+1 )+2*UP(K) )-F  PIE02550 

C  PIE02560 

C  *•»  FIND  RADIUS  FOR  (K+OST  APPARENT  WAKE  PIE02570 

C  PIE02580 

IF  ( (C«C-4»B*0) -LT.O. )  GOTO  80  PIE02590 

40  RP(K+1)»(-C+S0RT(C*«2-4-B-0) )/2./B  PIE02600 

.  IF(ABS(RP(N)-RT).LT. ( .0OO1»RT) )G0  TO  60  PIE02610 

IUP0AT=O  PIE02620 

C  PIE02630 

C  ••«  I0W=O  ->    TUNNEL  CORRECTIONS  PIE02640 

C         *\     =>  OPEN. WATER  (I.E..  INFINITE  FLUID)  PIE02650 

C  PIE02660 

IF(IOW.EO.I)  GO  TO  60  PIE02670 

IF  (ITER. NE- 1)00  TO  50  PIE02680 

ITER-2  PIE02690 

UENOLD=UE(N)  PIE0270O 

RPNOLD=RP(N)  PIE02710 

UE(N)=UE(N)-.98  PIE02720 

GO  TO  20  PIE02730 

50  CONTINUE  PIE02740 

ITER=ITER+1  PIE02750 

IF( ITER.GT. 10)G0  TO  80  PIE02760 

UEN=UE(N)  PIE02770 
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c 

c  •*• 

c 

60 


c 

c  •** 

c 


70 


C 

C  *■ 

c 


100 


c 

c  •** 
c 

110 

c 

c  •*• 

c 
c 

c 


UE(N)=U£(N)+(RT-RP(N) )  / ( RP ( N ) -RPN0LD ) * ( UE ( N ) -UENOLD ) 

RPNOLD=RP(N) 

UEN0LD=UEN 

GO  TO  20 

ITERATE  ON  HUANG  CORRECTION  FACTORS.  3  TIMES  FOR  OPEN  WATER 

IF( IUPDAT.EO. 1 )G0  TO  70 

IF( IOW.EO. 1 . AND. ITERAT.E0.5)G0  TO  70 

ITER=1 

ITERAT=ITERAT+1 

IF( ITERAT.GT. 10)G0  TO  80 

IUPOAT-1 

FIND  INOUCED  VELOCITIES  AT  APPARENT  WAKE  RADII 


CALL  EVALDK(NA.N.RA. RP.UAR.AUA) 

GO  TO  20 

CONTINUE 

CALCULATE  EFFECTIVE  WAKE  RADII 

NXM1=NX-1 

DO  100  M=M  .NX 

RE(M)=RH+( RPROP -RH)*FLOAT(M- 1 )/FL0AT(NXM1 ) 

IF  OUTERMOST  APPARENT  RADIUS  IS  INSIDE  PROPELLER  DISK.  THEN 
LINEARLY  EXTRAPOLATE  TO  PROPELLER  RADIUS 

IF  (RP(N).GE.RPROP)  GOTO  110 

UE(N)=UE(N1)  +  (UE(N)-UE(NI ) ) « ( RPROP - R ( N1 ) ) / ( RP ( N ) -RP ( N 1 ) ) 

RP(N)=RPROP 

INTERPOLATE  EFFECTIVE  WAKE  VELOCITIES  AT  SPECIFIED  RADII 

CALL  UGLYDK(N. 1 . 1 . RP . UE . 0. . 0. . AUE ) 
CALL  EVALDK(N.NX.RP.RE.UERX.AUE) 

COMPUTE  VOL  =  AVERAGE  VELOCITY 

V0L=V0L/3./(RT**2-RH»*2) 

COMPUTE  VOLUMETRIC  AVERAGE  NOMINAL  VELOCITY 

CALL  INTEOK ( NX , RX . RH . RPROP . YOX . XYDX . XXYDX . AUX ) 
V0LAVN=XYDX«2/( RPROP *«2-RH*«2) 


•-•  COMPUTE  VOLUMETRIC  AVERAGE  EFFECTIVE  VELOCITY 

CALL  I NT EDK(N.RP.RH, RPROP. YOX. XYDX. XXYDX. AUE) 

VOLAVE=XYOX*2./(RPROP-«2-RH-«2) 

RETURN 


PIE02780 

PIE02790 

PIE02800 

PIE028  10 

PIE02820 

PIE02830 

PIE02840 

PIE02850 

PIE02860 

PIE02870 

PIE02880 

PIE0289O 

PIE02900 

PIE02910 

PIE02920 

PIE02930 

PIE02940 

PIE02950 

PIE02960 

PIE02970 

PIE02980 

PIE02990 

PIE03000 

PIE03010* 

PIE03020 

PIE03030 

PIE03040 

PIE03050 

PIE03060 

PIE03070 

PIE03080 

PIE03090 

PIE03100 

PIE031 10 

PIE03120 

PIE03130 

PIE03140 

PIE03  150 

PIE03160 

PIE03170 

PIE03180 

PIE03  190 

PIE03200 

PIE03210 

PIE03220 

PIE03230 

PIE03240 

PIE03250 

PIE03260 

PIEO3270 

PIE03280 

PIE03290 

PIE03300 
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C  ...  TUNNEL  CORRECTION  OR  HUANG  CORRECTION  ITERATIONS  FAIL  TO  CONVERGE  PIE03310 

C  PIE03320 

80   STOP  PIE03330 

END  PIE03340 
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